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Abstract: We investigate the transverse- momentum-dependence in semi-inclusive deep 
inelastic leptoproduction of hadrons. There are two different theoretical approaches to 
study this dependence, one for low and one for high transverse momentum of the ob- 
served hadron. We systematically investigate their connection, paying special attention 
to azimuthal distributions and to polarization dependence. In the region of intermediate 
transverse momentum, where both approaches are applicable, we find that their results 
match for certain observables but not for others. Interpolating expressions are discussed 
for the case where one has no matching. We then use power counting to determine which 
mechanism is dominant in various azimuthal and spin asymmetries that are integrated over 
the transverse momentum. Our findings have consequences for the extension of transverse- 
momentum-dependent factorization beyond leading twist. They also shed light on the 
problem of resumming logarithms of transverse momentum for azimuthal distributions. 
Our results can be carried over to the Drell-Yan process and to two-hadron production in 
e + e~ annihilation. 
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1. Introduction 



Short-distance factorization is a key concept in quantum chromodynamics, providing much 
of the predictive power of the theory in high-energy scattering processes. Among the sim- 
plest processes to which this concept can be applied are one-particle or single jet inclusive 
production in lepton-nucleon scattering, two-particle or dijet inclusive production in e + e~ 
annihilation, and Drell-Yan lepton pair production via a photon or electroweak gauge bo- 
son in hadron-hadron collisions. Crossing of the hard-scattering subprocess closely relates 
these reactions, and many results obtained for one of them carry over to the other ones. A 
number of nontrivial issues for factorization arise especially when one observes the trans- 
verse momentum qr and the angular distribution of the produced particle with respect to 
a suitable reference direction. The problem then involves three scales, namely the scale of 
nonperturbative QCD dynamics, which we represent by the nucleon mass M, the trans- 
verse momentum qx, and the photon or electroweak boson virtuality Q, which throughout 
this paper we require to be large compared with M. 

There are two basic descriptions for the production of a particle with specified trans- 
verse momentum. One of them is applicable for qx <C Q and involves transverse-momentum- 
dependent (also called unintegrated) parton distribution and fragmentation functions. The 
other one requires that qx S> M and generates transverse momentum in the final state by 
perturbative radiation, using collinear (or integrated) distribution and fragmentation func- 
tions as nonperturbative input. In the following we refer to the two momentum regions 
and the associated theoretical descriptions as "Iow-c/t" and "high-g^", respectively. The 
low- and high-gr domains overlap for M <C qx <C Q, where both descriptions can hence be 
applied. An important question is whether in this intermediate qx region they describe the 
same dynamics or two competing mechanisms. Depending on the answer, one can either 
try to construct a formulation that smoothly interpolates between the two descriptions, or 
to add their results in a consistent manner. 

For the cross section depending on q\ but integrated over the angular distribution of the 
produced particle, the work of Collins, Soper and Sterman [1] showed that the descriptions 
based on intrinsic transverse momentum and on hard perturbative radiation indeed match 
at intermediate qx and permit a smooth interpolation at all orders in a s . A key element 
of the derivation was that for sufficiently large transverse momentum one can express 
unintegrated parton distributions and fragmentation functions in terms of their integrated 
counterparts and of perturbatively calculable hard-scattering kernels. The matching of 
the two descriptions allowed the authors of [1] to resum large logarithms ln(Q 2 /q T ) to all 
orders using renormalization group techniques — a procedure that remains the cornerstone 
for transverse-momentum resummation in a wide range of collider processes. 

For the angular distribution, however, the situation is less well understood. Both 
mechanisms just mentioned give rise to nontrivial angular dependence, as has been pointed 
out long ago for semi-inclusive deep inelastic scattering (SIDIS) by Cahn [2, 3] and by 
Georgi and Politzer [4] for the low- and high-g^ mechanism, respectively. To the best 
of our knowledge, the relation between the two descriptions for the angular distribution 
in unpolarized scattering has not been analyzed so far. Important progress has recently 
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been made in the understanding of a particular azimuthal asymmetry for a transversely 
polarized target in SIDIS or Drell-Yan production. The authors of [5-8] have shown that 
the description of this asymmetry by the Sivers effect for small qr and by the Qiu-Sterman 
mechanism for large qr match at order a s in the intermediate region M <C qr *C Q- It is 
natural to ask if one has a similar situation for other observables as well. 

In the present work we therefore present a systematic analysis of the interplay between 
the Iow-qt and the high-g^ mechanisms for angular distributions, both in unpolarized 
and in polarized scattering. This provides guidance for the theoretical description of a 
variety of observables, determining in particular whether or not one should add different 
contributions. For definiteness we will consider the case of SIDIS, but as remarked above, 
analogous studies can be performed for e + e~ collisions and for the Drell-Yan process. Our 
results are relevant to the possible extension of transverse-momentum resummation for 
specific azimuthal distributions, which was recently considered for the case of Drell-Yan 
production in [9] and [10]. 

A key finding of our work is that for certain observables, the leading terms of the 
low-and high-c/r descriptions match in the region M <C qx *C Q of intermediate transverse 
momenta, whereas for others they do not match. That this may happen can be understood 
already at the level of power counting. The Iow-c/t description, which uses transverse 
momentum dependent parton densities and fragmentation functions, is based on taking Q 2 
large compared with q T and all nonperturbative scales. We chose qr/Q rather than M/Q 
as parameter for power counting, since in the intermediate-^ region it is the larger of the 
two. Taking for example an observable F with mass dimension —2, we can thus expand 



where l n are dimensionless functions. In our applications, the term with index n will 
correspond to twist-n accuracy in the Iow-^t calculation, where n > 2. In the region of 
intermediate qr we can further expand the functions l n (M/qr) for small M/qr and then 
have 



with coefficients l n ,k- The high-gj- calculation, which is based on collinear factorization, 
treats both Q and qx as large compared with nonperturbative scales like M. The relevant 
parameter for power counting in the intermediate region is therefore M/qx, and we have 
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with dimensionless functions h n . In our applications, the term with index n will correspond 
to twist n in the high-gy calculation, where again n > 2. In the intermediate-gr region we 
can then expand h n (qT/Q) for small qr/Q '■ 



F(q T ,Q) = Jpl^K, 
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with coefficients h n ^. Since both ( |1.2| ) and ( |1.4| ) are valid in the intermediate region, we 
can identify the coefficients = hkn- As a consequence, a term of twist n in the low-gr 
calculation will only correspond to a term of the same twist in the high-g^ calculation if 
n = k. We will for instance encounter observables with 
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where the term with l 22 = ^2,2 is of leading twist in both the low- and high-g^ calculations. 
The term with l^ 2 = h 2 4 is subleading in the low-gr calculation and becomes subleading 
in the high-gy when one takes the additional limit qx <C Q. For other observables, we will 
find 
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(1.6) 



Here the term with ^,4 = ^4,2 is leading in the low-qx calculation but subleading in the 
high-^T one, whereas the reverse holds for the term with = ^2,4- The respective leading- 
order terms in the two calculations will hence not match in the intermediate region of qx- 
Which term in ( |1.6| ) is larger in given kinematics obviously depends on the relative size 
of the two small parameters qr/Q and M/qx- We will discuss in section [O] how one can 
construct interpolating expressions using both terms. 

An important question is which terms in the expansions ( |1.1| ) and (|1.3|), and hence in 
( |1.2| ) and (|P|), can actually be calculated in practice. We discuss this in some detail in the 
main body of the paper, but already mention here that in the high-t/T framework there is a 
large number of results at twist two and three. In low-q? framework factorization is rather 
well understood at twist-two level, whereas its status is less clear at twist three. Little is 
known about the validity of factorization at twist-four accuracy in either framework. In the 
example ( |1.6| ) one can thus envisage to compute the terms l 2 ,4 and /i2,4, which are leading 
in their respective power counting scheme. The simultaneous validity of the expansions 
Ql.2| ) and ( |1.4j) in the intermediate region requires that h^ 2 = l 2 ^ and l^ 2 = /i2,4, but 
at present one cannot check this explicitly because a calculation of the power-suppressed 
terms h^ 2 and l^ 2 is beyond the state of the art. 

Several investigations have been performed assuming factorization for twist-three ob- 
servables in the Iow-qt description. Detailed calculations at tree level [11-14] are found to 
be self-consistent and give results with a structure similar to that of twist-two observables. 
Their extension to higher orders in a s , including a proper treatment of soft gluon exchange 
has not been achieved yet, and the study [15] suggests that such an extension will not be 
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trivial. In sections |5.4| and we will investigate observables where the leading terms in 
the expansions ( |T^ ) and (|1.4| ) coincide and have the coefficient l^p = ^2,3. The twist-two 
quantity is readily computed, and its comparison with the result for 1%^ obtained with 
a candidate factorization formula will shed light on low-g^ factorization at twist three. 

In experimental analyses one often has to integrate over the observed qx in order 
to accumulate statistics. One may simply integrate an observable over q\ or consider 
weighted observables like f dq\ (q T /M) p F(q T ,Q) with some power p. If in turn the mea- 
surement of the ^T-dependence suffers from large uncertainties, then both a differential 
observable and weighted integrals will be affected with large errors, so that the simple 
integral J dq^ F(q T ,Q) may be the best quantity to consider from an experimental point 
of view. For the theoretical analysis it is important to identify the relative importance of 
the different qx regions in an integrated observable, and to clarify their interplay if several 
regions are important. Our results will allow us to address this question at the level of 
power counting. 

Our paper is organized as follows. In the next section we define the structure functions 
for SIDIS, which are the observables we study in detail in this work. To set the stage, we 
recall in section |3| some important results for SIDIS taken differential in qx but integrated 
over the angular distribution of the observed hadron, recalling in particular the foundations 
of qx resummation in this context. In section |I] we collect the well-known results of the 
calculation of SIDIS with qx S> M in collinear factorization at leading order in a s , and then 
approximate these results for qx <C Q. In section ||| we take the opposite path, recalling the 
results for SIDIS with qx <C Q and approximating them for qx ^> M. For this we need the 
behavior of distribution and fragmentation functions at high transverse momentum, and we 
will derive the corresponding power behavior of these functions based on general grounds. 
In section |6| we will see for which observables the calculations of the two previous sections 
match for intermediate qx and for which ones they do not. The consequences for integrated 
observables are discussed in section [/]. Whereas in section [| we derive the power behavior 
for all structure functions introduced in section ^, we give in section |8| explicit results 
for those observables that appear at twist two in the high-g-r regime. The comparison of 
the high-qr with the low-qx expressions will allow us to draw some conclusions about the 
unsolved problem of qx resummation for angular distributions, as well as the possibility of 
extending low-qx factorization to twist three. The main results of our work are summarized 
in section [9|, and some technical details are given in the appendices. 

2. Structure functions in semi-inclusive deep inelastic scattering 

The physical process we investigate in this work is semi-inclusive DIS, 

£(l)+p(P)^£(l') + h(P h )+X, (2.1) 

where t denotes the beam lepton, p the proton target, and h the observed hadron, with 
four-momenta given in parentheses. We allow for polarization of beam and target, but 
restrict ourselves to the case of an unpolarized final state, i.e. to the situation in which 
h has spin zero or where its polarization is not observed. The corresponding observables 
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Figure 1: Definition of azimuthal angles for semi-inclusive deep inelastic scattering in the target 
rest frame [23]. Phi and S±_ are the transverse parts of Ph and S with respect to the photon 
momentum. 



cover a variety of situations with different types of power behavior we wish to discuss. 
Many of them have been measured in experiment, see [16-21] and the recent review in [22]. 

Working in the one-photon exchange approximation, we define the photon momentum 
q = I — I' and its virtuality Q 2 = —q 2 . We use the conventional variables for SIDIS 

Q 2 P-Q P'Ph (00 , 

X= 2p—q> V = —V " = 1^7' M 

and write M and for the respective masses of the proton target and the produced 
hadron h. We take the limit of large Q 2 at fixed x, y, z, and throughout this paper we 
neglect corrections in the masses of the hadrons or the lepton. 

It is convenient to discuss the experimental observables for SIDIS in a frame where P 
and q collinear. We define the transverse part Pt, of P£ as orthogonal with respect to the 
momenta P and q. Likewise, we define the transverse part Sj_ of the spin vector of the 
target, as well as its longitudinal projection S\\ along P^. We further define the azimuthal 
angles (fth and (ps of P^ and with respect to the lepton plane in accordance with the 
Trento conventions [23], as shown in Fig. |l| Covariant expressions for the quantities just 
discussed can be found in [14] . Finally, we write A e for the longitudinal polarization of the 
incoming lepton, with A e = 1 corresponding to a purely right-handed beam. 

The lepton-hadron cross section can then be parameterized as [14] 

da a 2 y 

dx dy dz d<j)s d4>h dP 2 ± xQ 2 2 (1 — e) 

x {f UU ,T + eF VU ,L + y/2e(l+e) cos <\> h F%£+ h + ecos(2^) F™^ h 
+ A e y/2e(l-e) sin^F^ 



+ S\\ 



0! + sin 4> h Fffi* + e sin(2^) 



■sin 2<f> h 
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+ \S±\ 
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+ e sin(^ + <fe) i#* k+ * s) + e sin(3^ - <fe) F^ 1 



sin(30 h -0 s ) 



+ V2e(l+e) sm <t>s F^ s + ^(l+e) sin(2^ - fa) F s Jf^s) 



+ \S±\X £ 



e z cos 



S ) F%U">-M + V2e(l-e) cos 5 F™ S0S 



+ V2e(l-e) cos(2^ - F ^'^ 



LT 



(2.3) 



where a is the fine structure constant and e the ratio of longitudinal and transverse photon 
flux, 



l-y + y 2 /2- 



(2.4) 



The 18 structure functions F on the r.h.s. depend on x, Q 2 , z and P 2 ± and encode the 
strong-interaction dynamics of the hadronic subprocess 7* + p — » h + X. Their first and 
second subscript respectively specifies the polarization of the beam and the target. In the 
structure functions F UUT , F uu L and F^ T h ^ , Fy T ^ h the third subscript refers to 
the transverse and longitudinal polarization of the photon. 

To calculate the SIDIS structure functions it is convenient to use light-cone coordinates 
with respect to the directions of the relevant hadron momenta. We introduce light-like 
vectors n + and n_ with n + • n_ = 1 such that, up to mass corrections, n + is proportional 
to P and n_ proportional to P^. A rescaling 



n+ 



KTL A 



71- 



k n. 



(2.5) 



corresponds to boosts in the collinear direction. The off-collinearity of the process is de- 
termined by the vector 

q£ = (f + (l-r)xP fl - P£/z (2.6) 
with r = q^/Q 2 . For ease of notation we denote the length of this vector by 



q T 



-It Itu) 



1/2 



(2.7) 



so that q\ is positive. There is a simple relation between the transverse momentum qlt of 
the photon with respect to the hadrons and the transverse momentum P/\ of the produced 
hadron with respect to the photon and proton: Pt, = — zq T — IrzxP^ . The SIDIS cross 
section differential in q T instead of P 2 ± is hence equal to z 2 times the r.h.s. of ( |2.3[ ). 



3. Factorization and qx resummation 

In this section we recall some important results for the description of hard processes with 
measured qx, in particular the factorization for low qx formulated by Collins and Soper [24] 
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and its connection to the procedure of transverse momentum resummation by Collins, Soper 
and Sterman [1]. In the following we refer to these authors as CS and CSS, respectively. 
In the next two subsections we focus on the unpolarized SIDIS cross section differential in 
qT but integrated over the azimuthal angle (j)^. 

3.1 Collins-Soper factorization 

In the work of CS, factorization was derived for the production of back-to-back jets in 
electron-positron annihilation, or more specifically for e + e~ — * A + B + X , where A and B 
are two hadrons belonging to opposite-side jets in the e + e~ cm. In general the momenta 
Pa and Pb of the two hadrons are not exactly back-to-back because of their recoil against 
the additional particles X produced in the process. The cross section, or equivalently the 
hadron tensor, depends on the transverse momentum q T of the virtual photon w.r.t. the 
hadrons in the cm. of A and B, which is the analog of q T introduced for SIDIS in (]2.6|), 
For qx <C Q the CS paper derived a factorized expression of the hadron tensor, which is 
a convolution in transverse momentum of a soft factor U and two fragmentation functions 
D^/ a and Dfl a for the fragmentation of a quark or antiquark into A ot B. In addition 
there is a hard-scattering factor H, which does not depend on any transverse momentum. 
To be specific, Eq. (7.14) in [24] gives the following expression for the hadron tensor: 

cxTrl^^HF^^ci 72 ,^ 1 ^ 2 )! 2 E e ' / d 2 p T c?k T d% 

a J 

x 5^{p T + k T + l T - q T ) D^ a (z A ,p 2 T ; ( A ) D?/ a (z B , k T ; ( B ) U(l T ) + . . . , (3.1) 

where . . . stands for terms that either vanish after integration over the azimuthal angle of q T 
or are power suppressed in 1/Q. The index a runs over flavors of quarks and of antiquarks 
with fractional charge e a . For consistency within the present paper we have slightly changed 
notation compared with CS. 1 The individual factors H, D^ a , Df/ a , and U depend on an 
ultraviolet factorization scale fj,, which we have not displayed for brevity. The derivation 
by CS is done in an axial gauge specified by a spacelike vector v, and the dependence 
of individual factors on this vector is through the parameters Ca = — (2-Pa • v) 2 /v 2 and 
(b = — (2Pb • v) 2 /v 2 . More generally, (a and ( B serve as cut-offs for rapidity divergences 
and in a gauge-invariant definition of the fragmentation functions arise from path-ordered 
exponentials involving the vector v, see e.g. [25-27]. The hadron tensor as a whole is of 
course independent of Ca and Cb > so that the dependence on these parameters has to cancel 
between the fragmentation functions and the hard-scattering factor. The fact that the soft 
factor defined by CS does not depend on them is less obvious. 2 



1 We write Di instead of V for the fragmentation functions, v instead of n for the gauge fixing vector, 
and n+, n_ instead of va, «b for the light-cone directions. Our normalization condition for U differs from 
the one of CS by a factor (2tt) 2 . 

2 The four-vectors entering the construction of U are v, n+, u_ and It- The gauge vector v used by 
CS has a zero transverse component, so that the only scalar products involving v are v ■ n + , v ■ n_ and 
v 2 — 2(v ■ n + )(v ■ Gauges related by scaling v — * nv are equivalent, which leaves only a possible 

dependence on v ■ n+ /v ■ n_ . This is however excluded by boost invariance, which requires that U must not 



change under the rescaling (2.5) 
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In this work we will assume that the factorization of Collins and Soper also holds for 
the hadron tensor in SIDIS at low q T . Such an expression, albeit with some differences, has 
been obtained by Ji, Ma and Yuan [26], and another relevant investigation has recently been 
made by Collins, Rogers and Stasto [27], which gives us confidence that our assumption 
can be justified. The analog of (|3.1|) then reads 



^Dis^Trl^^T^l^siDis^C 172 ,^ 1 ^ 2 )! 2 ^ 6 ' / d 2 PT d 2 k T d% 

x 8W(p T -k T + l T + q T )f?(x,p 2 T] ()DUz,k 2 T ;( h )U(l 2 T ) + ... , (3.2) 

where one of the fragmentation functions has been replaced by the distribution function 
f? for quarks or antiquarks in the target. In the following we refer to the factorization 
expressed in (|3.2| ) as "CS factorization". The result ( |3.2| ) gives rise to just one structure 
function, 

F UU)T = |H(sC 1/2 ^ -1 Ch /2 )| 2 J> e « / d 2 p T d 2 k T d 2 l T 

x 8^ (p T -k T + l T + q T ) ft{x,&; C) Df(z, 4; ( h ) U(l 2 T ) , (3.3) 

where we recall that H, fx, D\, and U depend on a renormalization scale ji. For brevity 
we omit the subscript "SIDIS" in H from now on. 

In the intermediate region M <C qr <C Q, one can go further since at least one 



of the momenta pt, It i n (3-3) is of order and hence large compared with the 
nonperturbative scale M . For large transverse momentum the soft factor U (1%) can be 
calculated order by order in a s , whereas f\{x,Pj,',Q an d Di(z,k T ;Ch) are respectively 
given as convolutions of perturbatively calculable kernels with the collinear distribution 
and fragmentation functions ff(x) and Df(z). We will follow this path in sections || 
and ||. The accuracy of this procedure is however limited: up to mass corrections one has 
(COi) 1 ^ 2 = x _1 zQ 2 , and the power counting in the CS derivation requires both £ and (h 
to be of order Q 2 , so that the perturbative expressions for f^(x,pj,;Q and Df(z,k 2 ^;(i % ) 
involve large logarithms ln(Q/gr). Let us sketch how these logarithms are resummed in 
the work of CS. The variation of D% with £/, is described by the Collins-Soper equation, 
which gives dDi/dln^h as a convolution in transverse momentum of an evolution kernel 
with D\. Analogous considerations apply to the distribution function fi(x,p^; (). For the 
Fourier transformed functions 

/?(*, b 2 - = J d 2 p T e ib ^ f?(x,p 2 T ; () , 
D a ,(z, b 2 ; ( h ) = J d 2 k T e lbkT D^z, k 2 ; ( h ) (3.4) 
one obtains ordinary differential equations, whose solutions can be written as 

mx,b 2 ;C) = mx,b 2 ) exp[-£'(*C 1/2 ,&)] , 
DUz,b 2 ;C h ) = D1(z,b 2 ) ^ V [-S'{z-\ x h / 2 ,b)] , (3.5) 
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where the Sudakov factor S is constructed from the evolution kernel. The structure function 
in ( |3,3| ) can then be rewritten as 

Fuu,t= \H(Q,Q-fi)\ 2 J ^e-*^ exp[-2S'(Q, &)] U(b 2 ;^) 

xfJ(x,b 2 ;^Dt(z,b 2 ;fi) (3.6) 

with 

JJ(b 2 ;fi)= I d 2 l T e iblT U(%;n). (3.7) 

Here we have fixed the gauge parameters as x 2 ( = z~ 2 C,h = Q 2 and restored the dependence 
on the renormalization scale \x. The Sudakov factor S' resums large logarithms of Qb, which 
corresponds to large logarithms oiQ/qx in Fuu,T since the typical values of b in the integral 
Q3.6] ) are of order 1/qx- 

For b <C 1/M the factors S', U, /f and Df can be expanded in perturbation theory. 
To avoid large logarithms of fib in this expansion, one should take the renormalization scale 
of order 1/6. A common choice in the MS scheme is [i = bo/b with bo = 2e~~ <E ~ 1.1, where 
7e is the Euler constant. This simplifies a number of perturbative coefficients: in particular 
the 0(a s ) term in the soft factor is then zero, and one has U(b 2 , fi = bo/b) = 1 + 0(a 2 ) up 
to power corrections in Mb. The small-6 expansion for the distribution and fragmentation 
functions reads 

A a (x,6 2 ;/ u = bo/b) = Y^[d% ® fi)(x;n = h/b) , 

i 

z 2 Dl(z, b 2 ; [i = bo/b) = ^2(D{ ® Cjf) (z; f = bo/b) , (3.8) 

3 

where the indices i and j run over quarks, antiquarks and the gluon. f{(x; fj,) and D\(z; fj,) 
are the usual collinear distribution and fragmentation functions, and ® denotes the familiar 
convolution in longitudinal momentum fractions, 

1 dx 



{D®C)(z;/j)= C^D(-^n)c{z^). (3.9) 

./ <y Z \ Z ' 



With the scale choice \i = bo/b we find large logarithms of Qb in \H\ 2 = 1 + 0(a s ). These 
can readily be resummed using the renormalization group equation for this factor, which 
allows one to write \H(Q, Q; fj, = b /b)\ 2 = \H(Q, Q; fj, = Q)\ 2 e~^ Q ' b h In the intermediate 
region M <C qr Q one therefore has 

F UU:T =\H{^ = Q)\ 2 ^Y. xe l J (02 e " ib ' 9T exp[-5(g,6)] Ufa = bo/b) 



x 

i 
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with S = 2S' + R. Here we have used that for dimensional reasons H(Q, Q; fj,) depends 
on Q only in the combination Q/fJ,, so that the only Q-dependence in H(Q,Q;fj, = Q) is 
through the argument of the running coupling. An analogous statement holds for the b- 
dependence in U(b 2 ; fi = bo/b) at small b. The result ( 3.10| ) only involves the usual collinear 
distribution and fragmentation functions, together with factors H, S, U, C m , C out whose 
perturbative expansions are free of large logarithms. 



3.2 Collins-Soper-Sterman resummation 

We now turn to the region of large qx ^> M, where one can evaluate the hadron tensor in 
standard collinear factorization. To leading order in a s we have 



Fuu,T 



1 



Q 2 (2irz)< 



1 dx 
x 



1 dz (q\ (l-x)(l-z) 
z \Q 2 xz, 



it (?) «f (f) «r»> +/? © <§) ^r* 4) + n (?) * 



(3.11) 



with power corrections in M/qx- The hard-scattering coefficients Cuu,T for the indicated 
partonic subprocesses are functions of x, z, and q T /Q, and will be given in section |||. 



Approximating (3.11) for qx <C Q one obtains 



W = i ^2 E k(^) ) + fi(x) {Di ® P qq + D{® P gq ) (z) 



+ {P qq ®fl + P qg ®f 9 1 ){x)D a 1 {z) 



(3.12) 



with 



power corrections in qx/Q and in M/qx- The factor L is defined 



as 



L[ 9pj = 2C F \n%-?,C F 



q\ 



(3.13) 



and P qq , Pgq, Pgg are the DGLAP splitting functions at lowest order in a s , given in ( 4.2S| ) 
below. We see that a large logarithm of Q 2 jq\ appears in the fixed-order calculation 
when qx <C Q. Corresponding logarithms at higher orders in a s spoil the convergence of 
the perturbative series. Collins, Soper and Sterman [1] have shown that these logarithms 
exponentiate and that their resummation results in a factorized expression, which we will 
refer to as "CSS factorization". The discussion in the CSS paper is given for the cross 
section of Drell-Yan production differential in q T but integrated over the azimuthal angle 
of q T . The corresponding result for SIDIS is given in [28] and can be written as 



Puu,t 4E se «/ -^^'^ eM- S (Q> b )] 

x E ( C « ® fi) ^ = bo/b) £ {D{ ® qT) (z; \i = bo/b) . 



(3.14) 
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This form is valid for M <C qr *C Q- It can be extended to the full large-g-r region, qT S> M, 
by adding the difference between ( 3.1 1| ) and its approximated form (|3.12| ). For further 
discussion of this matching of resummed and fixed-order terms we refer to [29]. As an aside, 
we remark that at qT ~ Q the longitudinal structure function Fjju^l is parametrically of 
the same order as i*W,T> whereas at qx <C Q it is suppressed by a relative factor q T /Q 2 . 
One may hence also apply the CSS prescription to Fuu,T + £Fuu,L or to Fuu,T + Fuu,L 
instead of Fuu,T- The term to which resummation is applied is the same in all cases, and 
only the unresummed part of the fixed-order calculation is different. 

Let us see how ( 3.14| ) reduces to ( 3.12 ) at leading order in a s . The Sudakov factor 
S(Q, b) reads 



with 



where A\ 



S(Q,b) 



A (a.) = J2 A * 



Q da 2 T O 2 

1 A{a s ^))ln^ + B{a s ^)) 



b^/b 2 M 

7T 



(3.15) 



B(a s ) = Y, B k 



IT 



(3.16) 



fc=i fc=i 
Cp and B\ = —3Cf/2. The coefficient functions C m can be written as 

k 



C™{x;n = b /b) = 6*8(1 - x) + £C£W(s) 

fe=i 



in(fc)(V| (£[£ 
7T 



(3.17) 



and an analogous expansion holds for C out . Using the DGLAP equation we can evolve /i 
from the scale \i = bo/b to \x = Q and obtain 



/ 1 a (x;6o/^) = /i a (x;Q) 



27 



(P g9 «> /f + Pgg ® /f ) (*) In I® + O(o*) • 

°0 



(3.18) 



Evolving Z?i in the same way and putting everything together, we obtain 



Fuu,T 



d 2 b 



-ib-q T 



l-^C F (\n 2 b2Q 
2tt 



bH? 
ft 2 



3 In 



6 2 g 2 



mx; Q)-^{p m ® n + n« ® so (*) ln ^ + -f E (1) ® /?) (*) 

Dt{z; Q) - £ (-D" ® P M + Df P w ) (*) m ^ + ^ £ (i^ ® «) ( 



+ 0(a 2 ). 



(3.19) 



The running of a s is irrelevant at the accuracy of this expression. Expanding the square 
brackets one obtains a term /f (x; Q) Df(z; Q) of order a®, which is independent of b and 
hence gives a contribution proportional to 5^(q T ) to Fuu t x- Since we require qr S> M, 
this term must be discarded. For the same reason, the first-order coefficients and 
C^a do not contribute to Fjju,t at order a s . With the integrals [30] 



7 fro # ?t ■/ b o It 



(3.20) 
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we recover the lowest-order result ( 3.12 ) from ( |3.19| ). Going to higher orders in a s , 



one finds that the term with A\ in the Sudakov factor produces the leading logarithms 
a k m 2fc-l fQ2jq2 \ j n Fuu,T, whereas the next-to-leading logarithms a* ln 2fc ~ 2 (Q 2 /q^) also 
receive contributions from the coefficient B\, from the one- loop running of a s , and from 
the leading-order evolution of fi(x) and D\(z). 

The a s expansion of the CSS factorization formula Q3.14 ), which we have just performed 



to leading order, allows one to determine the functions S, C m and C ont at a given order in 
perturbation theory by comparing with the collinear fixed-order calculation in the high-gy 
region. The functional form of (3.14) was however derived by CSS using the result of CS 



factorization in the intermediate region M <C qx <C Q- This is not immediately obvious 
by comparing ( |3,14| ) with ( 3.10| ), because in the former expression there is no hard and 



no soft factor. As pointed out in [31], one can however introduce the hard factor into the 
CSS expression. \H(fi = Q)\ 2 depends on Q only through the argument of a s , so that 
one can use the renormalization group equation for the running coupling to rewrite it as 
\H(n = Q)\ 2 = \H{^ = b /b)\ 2 e R ( Q < b \ Since C in and C out in (Og ) are also evaluated at 



fi = b /b, we can combine factors into C" in = {H^ 1 C in and C' out = {H^ 1 C out at that 
scale, which gives 

F UUtT = \H(v = Q)\ 2 ±Y, xe2 a I (02 e " iMT e*P[S'(Q,b)] 

x E( C - in ® fi)fot* = b °/ b ) E( jD i ® C^ t )(z; f x = b /b) . (3.21) 

» 3 

with S' = S + R. Identifying C /in = U 1 / 2 C in , C /out = U 1 ' 2 C ou \ and S' = S, we finally 
recognize the CS result ( [3.10 ). 



The upshot of this discussion is essential in our context: the CSS derivation of trans- 
verse momentum resummation for Fjju^t and its analogs in e + e~ annihilation and Drell-Yan 
production makes use of two facts: 

1. CS factorization is valid for these observables at qx Q, and 

2. in the intermediate region M <C qr *C Q its results match those obtained from 
collinear factorization in the high-f^T region. 

In such a situation one can go further and construct expressions that interpolate between 
the factorization formulae for low and for high qx and are thus valid for all qx, from qT = 
to <7t ~ Q- A prescription for this had already been given by CSS, and a number of 
different methods have been proposed later; see [32] for a discussion and references. We 
shall not dwell on this issue here. 

3.3 Azimuthal dependence and polarization 

So far we have discussed only F UUT . The situation for the other unpolarized structure 
functions, F vu L , F^^ h , F^ 2 ^ h , cannot readily be inferred from the results of CSS. In 
section || we will see that the splitting functions appearing in the analogs of (3.12) for 
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^ SIDIS 



these three structure functions differ from the usual ones. Since the splitting functions are 
relevant at next-to-leading logarithmic accuracy, it is not clear if and how resummation 
beyond the leading logarithmic approximation can be performed in this case. An analogous 
observation for the angular distribution in Drell-Yan production has been made in [9]. 
(Resummation at leading logarithmic accuracy has recently been considered in [10].) 

The extension of CSS factorization for polarized scattering is relatively straightforward 
as long as one integrates over the azimuthal angle of Qj>. For SIDIS this concerns the struc- 
ture function Fll, and a corresponding calculation in this framework has been presented 
in [32]. CSS resummation for Drell-Yan production with longitudinal beam polarization 
has been investigated in [33,34]. A detailed discussion of polarization in the context of 
collinear factorization can be found in [35]. 

In the present work we will apply CS factorization to polarized scattering and to the 
SIDIS cross section depending on the azimuth 0^, as has been done in [36]. The factors 
tffi(x,pj<',C) and ^u^l O 2 ! ^Ti C/i) i n the factorization formula Q3.2j ) are then replaced 
by the quark-quark correlator $> a (x,p T ;() and the fragmentation correlator A a (z,k T ;(f l ), 
which will be defined in section |5.1| . The result reads 

oc \H(x( l/2 ,z- l (l /2 )\ 2 Y, e l [ d2 PT d 2 k T d 2 l T <5 (2) (p T - k T + l T + q T ) 

x Tr{$ a (x,p T ;()^A a (z,k T ;( h )Y} U(l 2 ) (3.22) 

and gives rise to a number of spin and azimuthal asymmetries at leading order in 1/Q. 
It is an open question if and how CS factorization can be extended to power suppressed 
observables, at least to those coming with one factor of 1/Q. The calculations of our work 
are relevant to this question, as we shall see in section jO| . 

A structure function that has received much attention in the literature is F^^p 1 
which arises when the initial hadron is transversely polarized. At low qx this observable 
is nonzero due to the Sivers effect [37]: the CS factorization formula ( |3.22[ ) gives a leading 
contribution in 1/Q proportional to the Sivers function f^p{x,p?r,) [12]. At high one 
can describe the same observable in terms of the Qiu-Sterman mechanism [38]. The cor- 
responding calculation uses collinear factorization at twist-three level, i.e., i^y^ is 
suppressed by 1/qr compared with Fjju,t- Calculating the behavior of the Sivers function 
at high transverse momentum, the analysis in [7, 8] has shown that at order a s the two 
descriptions exactly match in the intermediate region M <C qr Q- The situation is the 
same for the corresponding asymmetry in Drell-Yan production [5,6,8]. This suggests that 
for F^pp* ^ and its Drell-Yan analog it should be possible to use the CSS resummation 
procedure for large logarithms of Q/qT, as discussed in [39]. 

In the following sections we will derive the power behavior for the full set of SIDIS struc- 
ture functions, both in the low-gy and in the high-g^ description. This will in particular 
determine whether or not one can envisage to use CSS resummation for these observables. 
To determine the power behavior we can restrict our calculations to the leading order in 
a s . Since we will not attempt to actually perform a resummation of large logarithms, 
we need not go to 6-space as in (|3.10| ) or (|3.14 ). We will instead directly work with the 



momentum-space version ( 3.22j ) of CS factorization. In particular, we shall recover the 
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fixed-order result ( 3.12 ) for the intermediate region when expanding the CS expression 
|D| ) of F UUjT in the limit q T > M. 



4. Prom high to intermediate qr 



In this section we present the calculation of SIDIS structure functions at high transverse 
momentum in terms of collinear distribution and fragmentation functions, and then take 
the limit qx -C Q- We give explicit results for the six structure functions appearing at 
leading twist and order a s . For seven of the remaining structure functions, which are of 
higher twist or higher order in a s , there exist studies in the literature, which we will briefly 
discuss. 

In the high-gT calculation, the generation of large transverse momentum is described 
by hard-scattering processes at parton level. The diagrams for the contributions at first 
order in a s are shown in Fig. ^. We introduce the scaling variables 

Q 2 



X 



Pa 'Pb 
Pa ■ q 



(4.1) 



for the partonic subprocess, where p a is the momentum of the incoming parton and pb the 
momentum of the parton which fragments into the observed hadron h. We furthermore 
use the transverse momentum q^ introduced in (2.6). Neglecting mass corrections we have 
xp a = xP and pb/z = Ph/z, so that 



q£ = + (1 — r)xp^ —p%/z 
with r = q^/Q 2 . The partonic Mandelstam variables are then given by 



s = (q + Pa) 2 
u= {p a - Pb) 2 



l -^Q 2 



Q 2 ■ 



t=(q-Pbf 



1 - z 



Q 2 



1 — X 



(4.2) 



(4.3) 



Pa 




Ph 




(a') 



Pa 




(c) 



Pb 




Pa ~ Pb 

(b) 





Pa " Pb 



Figure 2: Feynman diagrams for the processes j*q — > qg, j*q — > gq, and j*g — > qq. 
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The structure functions defined in section g can be written as convolutions of hard- 
scattering coefficients with collinear parton distribution and fragmentation functions, 



„ 1 a s 9 f 1 dx f 1 dz ~( Qt (1 — x)(l — z) 



Q 2 (2ttz) 2 a J x x J z z \Q 2 xz 

^ (I) D Ki) c ^ q9) + f ?(l) D Ki) c &y* + # (I) D Ki) c ^r*° 

(4.4! 



as we have already seen in section 3.2. We recall that a runs over flavors of quarks and 
of antiquarks. Analogous expressions with different kernels C give the structure functions 
Fjjjj L , Fyy^ h , and F^-™ 2 . At order a s (but not at higher order) one finds the relation 

Fuu,l = IFuu^ ■ (4-5) 

The structure functions F LL and F^?^ h for longitudinal target and beam polarization 
are also given by expressions analogous to (4.4), with different kernels C and with the 
unpolarized parton densities /" and ff replaced by their polarized counterparts g\ and gf. 
The hard-scattering coefficients for the partonic processes *y*q — ► qg, ^y*q — ► gq, *y*g — ► qq 
can be computed from the respective diagrams (a, a'), (b,b'), (c, d) in Fig. ^, and those 
for 7*g — > gg, 7*5 — ► gg, 7*5 — > t/g are identical to their counterparts obtained by charge 
conjugation. Process by process we have 



7 q^qg 



( l + x 2 z 2 2 \ 

C UV , T = 2C F (1 - x)(l - z) + \ , (4.6) 

V xz q^J 

C%+» = -AC F [xz + (1 - x)(l - z)] ^ (4.7) 

qr 

C™* 2 ^ =4C F xz, (4.8) 

Cll = 2C F (2 (x + z) + 91) , (4.9) 

V xz q^J 



Q_ 

q T 



C co l^ = -4C F (x + z-l)^-, (4.10) 



1 q^ gq 



C UUiT = 2C F ( { l-Z)z + 1+ (4.11) 

\ z q F J 

C™^ = AC F [x (1 - I) + (1 - x) i] ^ ^ (4.12) 

z q F 

C™ 2<t>h = 4C F x(l-z), (4.13) 

C LL = 2C F ( 2x + 2(1 -z) + £2 + { \~ 2)2 , (4.14) 

V xz z q£ J 

C™^ =4C F (x-z)^^Q-, (4.15) 
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7 9 



C UUtT = 2T R [x 2 + (1 - xf] [z 2 + (1 - z 

1-x Q 



,2] 2 i£ 

! ^ g 2 



xz 



r ,cos 4>h 



cos 2<p h 
UU 



-4T R (2x - 1) (25 - 1) 



8T fi £(l-x) 



2 g T 



C ££ = 2T fl (2*-l) [^(l-f]y| 



£fCOS <f> h 

LL 



z q T 



(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 



with Cf = 4/3 and Tr = 1/2. The relation C uu L = 2C^ 2<^^, ' l holds for each individual 
subprocess. Our results agree with those in [32,40]. 

The behavior of the above results in the region qji, <C Q 2 can be obtained by rewriting 
the 5 function in Eq. (4.4) as [41] 



(q 2 (l-x)(l-z) \ 
\Q 2 xz J 



O 2 

6(1 - x)S(l -z)ln\ + 



6(1 -z) 



+ 



(1-5)+ 

where the plus-distribution is as usual defined by 

G(y) 



6(l-x) + 



illn^ , 
Q 2 q 2 T r 



dy 



^GM-Gjl) i 



1-1/ 



l-z 



(4.21) 



(4.22) 



We have written the hard-scattering coefficients in (4.6) to ( [1.20[ ) in a way that allows for 
an easy extraction of the leading power behavior at small qr/Q- The result is 



Fuu' 



L 01/3 \" 2 
2 27T 2 Z 2 2^ Xe * 



Qt 



ft(x) Dt(z)L(^)+ ft(x) (Dt ® P qq + D{® P gq ) (t 



+ {P qg ®fi+P qg ®f?)(x) Df( 



UU,L 



9 pcos 2<j> h 

lb uu > 



(4.23) 
(4.24) 



F, 



'COS 4> h 

UU 



1 



Qg T 2tt 2 z 2 



L V 9t / 

(4.25) 



7 cos 2<f> h 
UU 



1 a. 



+ (P> qq ®K + P' qg ®f 1 )(x)D<i(z) 
() 2 ,^ 2 E K « f/fWDf(^(?)+/fW {Dl®P^ + Dl®P'^)( 

+ (p; g ®/? + p^® 



(4.26) 
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and 



7 cos 4> h 
LL 



1 



Qg T 2n 2 z 2 



+ (AP gg ® 5 f + AP qg ® 5 ?) (x) 

J3 Frf(z) g a x {x) {pi 

L V <?t / 



(4.27) 



(4.28) 



The factor L contains a logarithm of Q 2 /?T an d ^ s gi ven m ( p-13 ) , and the convolutions 
are defined in ( |3.9|) . The splitting functions are given by 



c F 

C F 



1 + x 2 



_{l-x) + 
l + (l-5) 2 



T K [x 2 + (1 



f) 2 l 



2x 2 



(4.29) 



(l-x)+ 
2Ci7 (1 - z) 



^~gg (^0 



2T R x(2x-l) 



AT R x 2 



P>' q (z) = 2C F z 



for convolutions with unpolarized distribution or fragmentation functions, and by 
AP qq {x) = P qq (x) , ^Pqg(x) = T R (2x - 1) , 



AP'Jx) = PL 



qq\ 



^Pqg{x) 



2T R X 



(4.30) 
(4.31) 

(4.32) 
(4.33) 



for convolutions with polarized distributions. Similar results for the Drell-Yan process have 
been obtained by Boer and Vogelsang in Ref. [9]. We note that our P' qq and P qq correspond 
to P qq in Eq. (38) of Ref. [9] , while our P' qg corresponds to their P~ and our P qg to their 

P'qg/Z- 



cos 4> h 



Let us remark that the l/q T power behavior of Fyy 



and F 



cos 4> h 



LL 



arises as 



where \jq\ comes from the ^-channel propagators in the hard-scattering graphs of Fig. |2[ 
Likewise, the constant q T behavior (up to lng|, terms) of F UUL and F^ 2< ^ h arises as q\jq\, 
where the \jq\ from the hard propagators is fully canceled by numerator factors. We 
emphasize that, although they do not have a power-law divergence, the above expressions 
for Fjjjj L and F^ 2 ^ h cannot be used for small qx, because the approximations giving 
\jq\ for the hard propagators break down when gf<M. In fact, angular momentum 
conservation requires F^^ h and Fj^ to vanish like q T and F^ 2 ^ h to vanish like q\ for 
q T — ► 0, as shown e.g. in [42]. 



The six Fjjt structure functions for transverse target polarization vanish in the leading- 
twist approximation at high transverse momentum, because they would require the combi- 
nation of the transversity distribution function h\ with a chiral-odd collinear fragmentation 
function of twist two, which does not exist for an unpolarized hadron. They are however 
nonzero at twist three, where collinear quark-gluon-quark and three-gluon correlation func- 
tions appear, so that many more diagrams than the ones in Fig. Q need to be computed. 
Such a computation has been performed by Eguchi, Koike and Tanaka in Refs. [43,44]. 
The results for the different structure functions involve the product of D\ with Gf and Gf, 
which are chiral-even functions appearing in the decomposition of the quark-gluon-quark 
distribution correlator given in ( 5.37| ). Some observables involve in addition the product of 
hi with Ef, which is a chiral-odd function appearing in the decomposition of the quark- 
gluon-quark fragmentation correlator. The structure function Fy TT h ^ has also been 
computed in [7]. The result differs from the one in [43,44] because both calculations were 
missing certain terms. With the corrections discussed in [8], agreement between the two 
groups has been achieved. We note, however, that the twist-three calculation of the Fut 
structure functions is presently not complete. Terms involving Ep are only considered 
in [43], where they are found to contribute to E^ h ~* s) , F^ T (<ph+ * s \ and F§^ s . The 
calculation of that work is restricted to so-called derivative terms due to soft gluon pole 
contributions. The remaining soft gluon pole contributions, as well as contributions from 
soft fermion poles and hard poles are evaluated in [44], but only for Gf and Gp ■ In [8] 
it is shown that soft fermion pole contributions from further diagrams must be added to 
those results. Finally, all calculations in the literature are restricted to quark-gluon-quark 
functions of twist three, so that three-gluon correlators do not appear. 

Using the soft gluon pole and the hard pole contributions computed in [43,44] we have 
extracted the leading behavior of all Fjjt structure functions in the limit <C Q- The 
structure of the results is listed in Eqs. Q6.8| ) to ( |6.13|) of section ^, both for the power 
law and for the distribution and fragmentation functions appearing in each observable. 
We have verified that this structure is not changed by the soft fermion pole contributions 
given in [44]. Since the corrections to [44] discussed in [8] concern only soft fermion pole 
contributions in p^^ h+( ^ s ^ they do not affect the structure of (|6.10| ) either. The same 
should hold for the remaining five Fjjt structure functions, but this has not been checked 
explicitly. 

The structure function F^ is nonzero at twist two and order a 2 s . Depending on a 
single polarization, it is a T-odd observable and hence requires an absorptive part in the 
amplitude, which in this case is provided by a loop in the hard-scattering subprocess. The 
relevant graphs have been calculated in [45], and numerical estimates for specific kinematics 
have been given in [45,46]. For the structure functions F^^ h and F^ 2< ^ h the situation 
is similar, but no explicit calculation exists in the literature. There is no contribution to 
psm<p h ^ phincph ^ jpmn2<t> h ^ th ree anc [ orc i er because the necessary T-odd terms 

would need to come from twist-three quark-gluon-quark correlators. For an unpolarized or 
longitudinally polarized hadron, these are chiral-odd [43] and have no twist- two chiral-odd 
partners in the other correlator. From the calculation in [45] we can extract the power 
behavior of F^^ h for q T <C Q, which we will give in Eq. ( |6.5| ). 
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5. Prom low to intermediate g^: power counting 



In this section we derive the behavior of distribution and fragmentation functions at high 
transverse momentum. Plugging the results into the known Iow-^t expressions of the 
SIDIS structure functions, we will obtain their power behavior in the intermediate region 
M <C qr "C Q. To begin with, we specify in the next two subsections the distribution and 
fragmentation functions that will appear in our calculation. 

5.1 Transverse- momentum-dependent distribution and fragmentation functions 

For the discussion of distribution and fragmentation functions we use light-cone coordinates 
defined with respect to the momenta P and Ph, which we already introduced at the end 
of section |2[ For any four- vector a we then have the plus-component o + = a • the 
minus-component a~ = a ■ n+ , and the transverse part a!t = a M — a + n^_ — a^n^. The 
hadron momenta read 



M 2 „ „„ „_ „ Mj 

and the spin vector of the target can be decomposed into longitudinal and transverse 
components as 

P+ ^_M_ 
M U+ 2P^ 

The transverse-momentum-dependent quark distributions appearing in the description of 
SIDIS are defined from the quark-quark correlation function 



S» = S L (—n1-—nt)+S£. (5.2) 



(5.3) 



where p + = xP + and summation over color indices is understood. The corresponding 
correlation function for antiquarks is obtained by replacing the quark fields by their trans- 
forms under charge conjugation, see Ref. [11]. The quark fields in (5.3) are renormalized 



fields, and the corresponding scale dependence of the correlation function is given by a 
renormalization group equation involving the quark anomalous dimension [24]. 

The gauge link U(n^\ in (5^) is a Wilson line that connects the quark fields and thus 



makes the correlation function color gauge invariant. In the factorization theorems for 
scattering processes, the gauge link incorporates the exchange of gluons between partons 
that move in the opposite light-cone directions n + and n_. Consideration of gluons collinear 
to the target yields Wilson lines with paths that point along n_ and lead to light-cone 
infinity, a~ = ±00, where they are closed by transverse segments from Ot to £t [13,47]. 
Different processes require different gauge links. In particular, the simplest links closed at 
tz~ = ±00, which we denote by U^, give rise to the correlators & + \x,pt) and &~\x,px) 
appearing in SIDIS and Drell-Yan production, respectively. More complicated gauge links 
show up in other processes [48,49]. 

When defined with strictly lightlike Wilson lines, the correlator ( |5.3[) contains diver- 
gences in gluon rapidity (sometimes referred to as "endpoint singularities" ) and hence must 
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be modified [25]. Different schemes have been discussed in the literature. One possibility 
is to use paths that point in a non-lightlike direction v instead of n_ [26,27]. Up to subtle 
issues we will mention in Appendix |A], this is equivalent to working in axial gauge, A-v = 0, 
which was used in the original scheme of Collins and Soper [24]. In a number of different 
schemes, the proton matrix element in ([TI^) is divided by vacuum expectation values of 
suitably chosen Wilson lines [27,50-52]. The arguments in the present section use Lorentz 
invariance and power counting, so that we need not specify the detailed choice of scheme. 
As long as v is a linear superposition of n+ and n_, no new four- vector is introduced in 
pt), which therefore depends on v only via the scalar parameter £ = —(2P ■ v) 2 /v 2 we 



already encountered in section 3.1 



The correlation function ( |5.3|) can be parameterized in terms of distributions functions 
depending on x and p\ as [14] 

= \ j/i 4+ + fir Tp If 70 4+ + 9ilS l ^5 4+ ~ 9it ST jJ T 75 4+ 

. , 75 [#r, 4+] i j_ s t p PtPt 75 [>, 4+] , i I c 7s|>T, 4+] . u ±i[h, 4a 
+ n l 7, n lT 773 n r H lL b L — — h % 



M J . S T -p T . ± S Tp e p T p Ta 



M 2 2 ijL ^ 2M 1 2M 



2P+ I M M 

± u i f . j> a per 
+ f ~7T + /l 77 7ct ~~ /t TTo 7 + h Sr P ^ la 



pa o \P u 1 

T , ,l n PTp^T f ± ^TpPrPT e Tar T , Q pa 

■7 +JT b Tp e 

, l c 75j^T ±PTp e T ± S TpPTPT . a 

+ 9lS l ^~ + g M 757a - ffr ^2 757<t + 0t75#t 

± 75 [ *$T; jfej 5 T - Pr 75 [^+, 75 [^+, ^_] , gfj 

T 2M T M 2 L L 2 2 

M 2 r , 

where the two-dimensional antisymmetric tensor is given by 

= e al3pa n +p n- a (5-5) 

with e 0123 = 1. Index pairs in parentheses indicate that the trace is subtracted in the two 
transverse dimensions, 

PtPt = PtPt ~ o ( p T " Pt) 9t > ( 5 - 6 ) 
where the transverse metric tensor is g^P = g al3 — rP_ — n°_ . The first eight distribu- 



tions in (|5,4j ) are referred to as twist two, and the next sixteen distributions as twist three. 
The . . . stand for the remaining eight distributions of twist four, which are given in [53]. 
We will not need them in the following and tacitly omit them in further parameterizations. 
Corresponding to the Dirac matrix structure in the decomposition ( |5.4| ), functions denoted 
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Table 1: Behavior of distribution functions under time reversal and in the high-p^ limit. The 



time reversal factor rjf is defined in (5.7) and the exponent n for the high-p-r behavior in (5.8 



with letters /, g or e, h are respectively referred to as chiral-even or chiral-odd. Functions 
with subscripts L or T appear in the parts of $>(x,pr) that depend on the longitudinal 
or transverse component of the spin vector. (An exception to this rule of notation is the 
transversity distribution hi.) 

It is understood that the correlator and each of the functions in (|5.4D should carry 
a label specifying the gauge link, as well as a label for the quark flavor. Time reversal 
connects with & uT \ and in particular $W with <I>H. This provides relations [54] 

fW(x,p 2 T )= Vf fM(x,p 2 T ), (5.7) 



where / stands for any of the distributions in ( |5.4| ). We call a distribution T-even if rjf = +1 
and T-odd if r\f = —1. The values of rjf are given in table [5]. We also anticipate in the 
table the power behavior 

fW(x,p 2 T ) ~ (5.8) 



of the distributions for px 3> M, which we shall derive in section 5.5. 
Fragmentation functions are defined from the correlator 



A.(^t) = ^E/ 



d£ + d 2 £ T eik< 



(27T) 



3 



X 

out out 



(5.9) 



f-=0 



where k~ = PZ jz and N c = 3. The prefactor 1/(2N C ) comes from averaging over the 
polarization and color of the fragmenting quark. The subscript oo in the gauge links 
indicates a space-time point with plus-coordinate a + = oo. The precise choice of Wilson 
lines involves the same issues we mentioned for the distribution correlator. Aspects specific 
to the case of fragmentation functions are discussed in [49,50]. Notice that fragmentation 
functions with different gauge links are not related by time reversal, because time reversal 
transforms "out" states \h,X) out into "in" states \h, X){ n , and the difference between these 
states amounts in general to more than just a phase. 
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For an unpolarized hadron h the decomposition of the fragmentation correlator reads 



+ 



2Pr 



M h T 2 



k 



Tp £ T 



757o 



(5.10) 



where the functions on the r.h.s. depend on z and kj,. In a more explicit notation they 
should also carry a flavor index. 



5.2 Collinear distribution and fragmentation functions 

In applications of collinear factorization, the structure of incoming hadrons is represented 
by the light-cone distribution correlator 



2vr 



(p\^(o)u?- MO\P 



(5.11) 



=0,£ T =O T 



where the gauge link U^g, connects the quark fields along a path in the n_ direction. This 
would simply be the integral of the p^-dependent correlator introduced in the previous 
subsection, 

$(x) = [ d 2 p T <S>M(x,p T ), (5.12) 



were it not for two complications. On the one hand, the correlator ( 5.11] ) has ultraviolet 
divergences due to the fact that all field operators are taken at the same transverse position. 
Their regularization and subtraction gives rise to a scale dependence described by the 
DGLAP equations. Correspondingly, the integrand on the r.h.s. of (5.12) diverges like l/pj> 
at high p T , as we will see in the next subsection, so that the p T -integral must be regularized. 
On the other hand, the rapidity divergences of the pr-dependent correlator, which we 
discussed in the previous subsection, cancel under the integral over px [25,51]. They require 



no regularization in the collinear correlation function ( 5.11 ), which hence is independent 
of the parameter £. The different regularization procedures in the correlators $>(x,pr) 
and $(x) reflect the different types of subtractions required when constructing transverse- 
momentum-dependent or collinear factorization theorems. We will shortly discuss how the 
relation fl5.12| ) should be understood. 



The correlation function (5.11) can be parameterized as 



*W = 2 j/i i+ + 9iS L J$ 4+ + th > + | e - e L S L ij$ 

, f c P° , <r , u a 75[^+, 44 . ,i[4+,4-] 
+ /T S Tp e T 7 ff + 5t75 $t + h L S L + h 

where the distributions on the r.h.s. depend only on x. They are given by 



hi?) 



j d 2 p T /i(s,Pr) 



(5.13) 



(5.14) 
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and similarly for the other functions, with one common exception of notation, 

gx(x) = j d 2 p T g 1L {x,p 2 T ) . (5.15) 



Since the gauge link in (5.11) is along a finite light-like path from to £, time reversal 
relates the collinear correlator with itself, and as a consequence fr(x) = ei{x) = h(x) = 0. 
This ensures that ( |5,12| ) can simultaneously hold for different links U, in particular for U + 



and IA , which according to (5.7) give p^-dependent distributions and of 

opposite sign. 



To make the meaning of ( 5.12 ) more precise, we observe that the combination of 
Eqs. and (fy|) gives 



/l(x,6;C,M) = J d 2 p T e ibp Tf 1 (x,p 2 T -X,v) = fi(x-,fi) + 0{a s ln 2 ((b 2 )) (5.16) 

for small enough b, where we have set fj, = bo/b and used the perturbative expansions 
S' = 0(a s ln 2 ((b 2 )) and 6% = 5 ai 5{l - x) + 0(a s ). The factor e ibp T in ( p6| ) regulates 
the logarithmic divergence of / d 2 p^ fxi^iPr) by damping the integrand for large pt > \/b. 
Alternatively one can cut off the integral at pr = bo/b, since 

j d 2 p T e ib - p T / 1 (x,^;C,/i) = J d 2 p T 0{is 2 -p 2 T )h(x,p 2 T] C,ri (5.17) 

up to corrections of order b 2 , as we will show in appendix We thus see that the relation 
(5.12) should be understood with a suitable regulator of the integral on the r.h.s. and as 



up to corrections of order a s . The same holds for ( |5.14| ), ( |5.15 ), and for similar integral 



relations in the following. Let us remark that an extension of ( 5.16j ) to the full correlation 
function <I> has not been given in the literature. 

As we will see in the next subsection, calculations at subleading power or those in- 
volving azimuthal asymmetries lead in the collinear expansion to p^-weighted correlation 
functions 

< ±] (x) = J d 2 p T p^^(x, PT ), (5.18) 

where the Lorentz index a is restricted to be transverse, and where the same remarks 
about regularization apply as for ( |5.12 ). In contrast to &(x), the correlator $g^(a?) does 



depend on the choice of Wilson lines in &^(x,pt) and hence contains both T-even and 
T-odd distributions. Omitting the superscript [W] for the sake of legibility, we have the 
decomposition 



+hi p SL j ! bMA +h m'M^Ay^. { ... } , (5 . 19) 

where we have only displayed the terms of leading twist and defined p\ moments 

// 2 \ n 

d 2 PT{02) f&M), (5-20) 
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and similarly for the other functions. The functions f^f 1 and are T-odd and thus 

change sign when going from *gM(a?) to 

The factor pj- in (|5.18|) can be converted into a derivative d a acting on the matrix 
element that appears in the definition ( |5.3| ) of $>(x,pr)- One can then express <&g(x) 
in terms of correlators with either a gluon field or a covariant derivative between the 
antiquark and quark fields. The former is a collinear quark-antiquark-gluon correlation 
function, whereas the latter can be rewritten in terms of the quark-quark correlator <J?(x) 
using the equation of motion for the quark field. In this way, the moments given in 
J5.19 ) can be traded for functions of twist three, up to twist- two distributions multiplied 



by the quark mass [14]. 

The /c^-integrated fragmentation correlator for an unpolarized hadron has the decom- 
position 

A(z) = z 2 J d 2 k T A(z, k T ) = \Dii- + ^ {E + H l[i - 2 i+] } , (5.21) 

where the fragmentation functions on the r.h.s. depend on z. They are given by 

Dx(z) = z 2 J d 2 k T D±(z,kl), (5.22) 

and similarly for the other functions. Notice the factor z 2 , which appears because D\(z, k 2 -,) 
is a probability density w.r.t. the transverse momentum k' T = —zk T of the final-state hadron 
relative to the fragmenting quark [11,55]. As already discussed, time reversal invariance 
does not constrain fragmentation correlators, so that H(z) can be nonzero unlike its distri- 
bution counterpart h(x). For the fer-weighted correlation function needed in calculations 
at twist three and higher, we have 



ASW^'/^^AW^-^^'^ + lf...) (5.23) 



with 



Ht {n \z) = z 2 f d 2 k T (Jfe) (5-24) 



where again a is restricted to be transverse. 

5.3 Distribution and fragmentation functions at high transverse momentum 

We are now ready to derive the behavior of correlation functions at high transverse mo- 
mentum. We consider the distribution correlator $(x,pr) for transverse momentum pt 
much larger than the scale of nonperturbative interactions. The generation of the large 
transverse momentum can be described in perturbation theory. Technically, we approxi- 
mate $>(x,pr) in powers of 1/pr using a collinear expansion that leads to the factorization 
of the transverse momentum dependence. To derive a formal proof of factorization, one 
would use the same techniques as for, say, the production of a high-pr jet in deep in- 
elastic scattering. We shall not attempt this here, but limit ourselves to determining the 
power-law behavior of the distribution functions that parameterize $>(x,pr), using Lorentz 
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(a) (b) 

Figure 3: Example diagrams for the calculation of the high-p^ behavior of the quark-quark 
correlator $>(x,pr)- The dashed lines represent the final-state cut. 

invariance and dimensional analysis as our main tools. The explicit calculation in section || 
and the one for the Sivers function f^p in [6,8] provide examples for the consistency of the 
collinear factorization formalism at leading order in a s and to leading and first subleading 
power in 1/pt- One should, however, be aware that factorization might break down at 
some higher-order or higher-power accuracy. 

The evaluation of &(x,pt) at lowest order in 1/pt involves diagrams as in Fig. ^a, 
whereas at higher orders correlators with three or more partons appear at the bottom of 
the graphs, as shown in Fig. ||b. To set up the power counting, we generically write p for 
the hard scale and take 

p + ,p-,p T ~ P+~ l + JtJt ~ P, (5.25) 
ItMtMt ~ M, (5.26) 
P"~ r,l^,l 2 ~ M 2 /p, (5.27) 



where the mass M represents the soft scale. In ( 5.27 ) we have used that the loop momenta 
I, l\ and l 2 are attached to a soft function at the bottom of the graphs and thus have virtu- 
alities of order M 2 . Starting point of the calculation is the correlation function depending 
on the full four-momentum p, 

*f(v) = / (04 eiP< MO\P) , (5-28) 

from which we obtain <&^(x,pt) = J dp~ &^(p). We omit the superscript [±] for clarity 
of notation in the next few steps, and will restore it when required later on. We restrict 
ourselves to the leading and first subleading order of $>(x,pr) in the 1/p expansion. To this 
order, the relevant factorized graphs can be written as the convolution of hard-scattering 
kernels and correlation functions in the form 

*(x,Pt) = J dp-dH £r(p,0*(0 + J dp-dH 1 d 4 l 2 H%{p,l 1 ,l 2 )^ Aa (l 1 ,l 2 ) 

+ {terms with two- and three-gluon correlators} + . . . (5.29) 

with . . . representing terms of higher order in 1/p. It is understood that the hard-scattering 
kernels H(p,l) and H^(p,li,l 2 ) include 5 functions putting the cut lines on shell — this 
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can readily be used to perform the integration over p~ . The lower blob in Fig. ||b is 
parameterized by the quark-gluon-quark correlator <£^t, which contains the gluon potential 
A a between the quark and antiquark fields. The gluon polarization index a in ( |5.29 ) is 



restricted to be transverse: the contribution from A~ gluons is power suppressed by at 
least 1 jp 2 , whereas the corresponding contribution of A + gluons ends up in the gauge link 
of the quark-quark correlator when the terms in the factorization formula are arranged in 
a gauge-invariant manner. 

We now expand the hard-scattering kernels in the small momentum components ( 5.26| ) 



and ( |5.27| ). To the order we are considering, we can neglect , l 2 and ^\ r p^ 
H^ip, h, h), whereas in H(p,l) we can neglect l~ but must expand the /^-dependence 
up to first order. This gives 

/ dp-H(p,l) = ^-xH 2 (x, P +, PT ) + l -^xH«(x, P + ,p T ) + ... , (5.30) 
J P + P + 

J dp~ H%(p,h,h) = ^x 1 x 2 H% j3 (x 1 ,x 2 ,p + ,p T ) + ... , (5.31) 
where we have introduced the plus-momentum fractions 

x = p + /l+, Xl =p + /lf, x 2 =p + /l£ (5.32) 

and chosen prefactors such that H 2 , H 3 and H^ 3 are invariant under longitudinal boosts, 
i.e. under the rescaling ( |2.5| ) of n± and the corresponding change of plus- and minus com- 
ponents. The convolution Q5.29 ) now takes the form 



Hx,p T ) 

= J ^H 2 (x,p + , PT ) J dl-d 2 l T m + J ^H«(x,p + , PT ) J dl-d 2 l T lTaMl) 

/dx dx f* 
^ H%Ax 1 ,x 2 ,p + ,p T ) / dli dl 2 d 2 l 1T d 2 l 2T <f> Aa (h,l 2 ) 
x x x 2 J - 

+ {terms with two- and three-gluon correlators} + . . . . (5.33) 
In the first two terms we recognize the collinear quark-quark correlators from ( |5.12 ) and 



$(y) = J dr d 2 l T $(/) , $g(y) = J drd 2 l T lT$(l) (5-34) 

with y = l + /P + , whereas the third term involves a collinear quark-gluon-quark correlation 
function 

= / dl^dl 2 d\ T d 2 l 2T n(h,h) (5.35) 

with y l = if /P + and y 2 = l 2 /P + . In order for these correlators to be gauge invariant 
one must reshuffle certain pieces among the different terms in ( 5.33j ), as shown for instance 
in [13]. On the r.h.s. of ( |5,34j ) and ( |5,35| ) this implies subtraction of terms with the gluon 
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potential A a at light-cone infinity, which we have not displayed. One also finds that taking 
the gauge link in $>^(x,pr) leads to the corresponding path-dependent correlators 
$ a d [±] (y) and 

i $g(yi,y 2 ) 



iP+ yi-y 2 ± ie 



(5.36) 



in the factorization formula, where 



2ir 2vr 

x (P|^(0)W ( ^ 2) gG +a fa)lfe Ai) MZi)\P) 



(5.37) 



does not carry a superscript [±] because, like &(y), it involves only Wilson lines of finite 
length along n_. We now decompose the correlation functions into terms of definite twist, 



M 

§{y) = ®2{y) + p+ *a(y) + • • • , 

M 

$A(yi,V2) = -p+ *5,3(ift s w) + -- - > 



$g( 2 /) = M$g )3 (y) + ... 



(5.38) 



where the prefactors are chosen such that $3 and $9,3, &a,3 are dimensionless and in- 
dependent of P + . Under a longitudinal boost $3 is invariant, whereas the other correlators 
transform like n + . Dimensional counting readily gives H 2 ~ 1/p 2 and H^,Ha,3 ~ 1/.P 3 - 
Using p + = xP + and ( |5.32| ) we can then rewrite (5.33) as 

dx 

x 



/dx / x \ 

— H 2 (x,p + ,pt)<S>2[-) 
x \xJ 



+ M { I — 

X 



+ 



dx H 2 (x,p + ,p T ) 



,-<!>,(-) - / ^H^p+^^J^ 



H A,?,{Xi,X 2 ,p + } p T ) X^ A ,3a(^, — 

\Xi X 2 



/dx\ dx 2 
xi x 2 

+ {terms with two- and three-gluon correlators} + 0(l/p ) 



(5.39) 



where the first term is of order 1/p 2 and the terms with prefactor M are of order 1/p 3 . 
To obtain the high-pr behavior of the individual distribution functions parameterizing 
$>(x,pr), we need to analyze the dependence of the hard-scattering kernels on p + and pr- 
The kernels carry four Dirac indices, so that ( |5,39| ) explicitly reads 



/dx / x \ 

-T- H 2tijk i(x,p + ,p T ) $2,fcz(^J +••• , 

and similarly for the terms of order 1/p 3 . We can decompose H 2 as 



(5.40) 



H 2 (x,p + ,p T ) = -j 
Pt 



^2 T m ® r n t mn (x,p + ,p T ) + ^ r ^-M (g) V n,v ^ni^ , P + , Pt) 



+ a XfM a up t x ^(x,p + ,p T ) 



(5.41) 
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with T m G {1,75} and T m ^ £ {7^7^75}, where the first matrix in the tensor products 
carries Dirac indices ik and the second one indices jl. Only an even number of 7 matrices 
appears in this decomposition, i.e., no structures like T m ^ (g) Y n or T m ^\ <g) F n>fll/ , which 
reflects that chirality is conserved in the hard scattering kernel. The scalars t mn and the 
tensors t% n , t x ^ p are dimensionless and invariant under longitudinal boosts, and therefore 
they can be constructed from g^ u , e\^ up , and the vectors 

U 4. U I I u 

Ur p^n\_ \p T \n , , . 

IKTI . (5.42) 



\Pt\ \Pt\ ' P + 

Since the tensors have an even number of indices, the factors \pr\ combine such that 
H2(x,p + ,pt) depends only on integer powers of \jp\- The same is readily shown for 
the kernels H^(x,p + ,pr) and H\ 3 (x 1 , x 2 ,p + ,Pt), which go like l/p 3 instead of l/p 2 but 
involve one more Lorentz index in the analog of the decomposition (5.41). Analogous 
arguments can be given for the kernels connected with two- or three-gluon correlation 
functions at the bottom of the graphs (which have two instead of four Dirac indices and 
two additional Lorentz indices for the exchanged gluons) and for the kernels that appear 
when &(x,pt) is evaluated to order l/p 4 or higher. 

The upshot of this argument is that the distributions parameterizing the correlator 
Q(x,pt) behave like integer powers of l/p\ for pj, M 2 . Together with the constraints 
from dimensional counting and Lorentz invariance, this allows us to determine the leading 
power behavior of each distribution. Matching the dependence on p + , we find for instance 
that terms on the r.h.s. of ( |5.39| ) contribute to the twist-two and twist-three parts of 
*(x,pr) = $2(37 Pt) + (M/P + ) ^> 3 (x,p T ) + ... as 

1 p p 

Pt P + Pt 

M , . Mp p T . . t , 

— -y -» x$ 3 (x,p T ), -» $ 2 (x,pr), (5.43) 

where we have used that 1/P + = x/p + . Comparing with the parameterization ( |5.4| ) of 
Q>(x,pr) we see e.g. that f£r(x,p T ) and gix{x,Pj<) decrease as M 2 /p T . If there had been 
terms going like p T /\p T \ 3 on the r.h.s. of ( |5.39|) , they would instead decrease as M/\p T \ 3 . 

A number of selection rules specify which collinear distributions can contribute to 
the high-pr behavior of a given p^-dependent distribution. Clearly, the dependence on 
the target polarization must match. Because the hard scattering conserves chirality, the 
chirality of distributions must match as well. Finally, we recall that the correlator <&{x,pr) 
depends on the gauge link and contains terms which are even or odd under the exchange 
<J>M(x,£>t) <-> <I>H(x,pt)- The T-odd terms in <&^\y) and $^(2/1,2/2) thus contribute 
to the T-odd distributions in &^(x,pt), and vice versa. The collinear correlator &(y) 
only contains T-even terms, but it can contribute to T-odd distributions in $W(x,px) 
through graphs with absorptive parts in the hard-scattering subprocess, starting at order 
a 2 . An example is shown in Fig. |j. Graphs involving the quark-gluon-quark correlator 
(such as the one in Fig. ||b) have absorptive parts already at order a s , which provides 
further contributions to the T-odd part of ^^(x,pr)- 
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Figure 4: Example of a diagram at order whose absorptive part gives rise to T-odd terms in 
the quark-quark correlator. The double line denotes an eikonal line originating from the gauge link 
in §>(x,pt), as specified in appendix |a|. 



With the time reversal properties given in table [j], we find that at order a s collinear 
functions denoted by letters /, g, e, h can only contribute to the high-pr behavior of 
distributions denoted by the same letter. Putting everything together we find 

fi ~ -5- asFlfi) , g lL ~ -o-a s Flgi] , hi ~ -5- a s F [hi] , 

p T Pt Pt 

xf L ~ \asF\fi] , xgj; ~ \ a s T[gi] , 

Pt P t 

xh T ~ — g- a s J 7 \hi\ , xJit ~ — a s F\h\\ (5.44) 

Pt Pt 

from the 1/p 2 part of $>(x,pr), and 

/-L M ' TtfMi) 1 m2 <r\ (i) 1 

J it r r a sF[fiT >•••]) Sir r r a « F[9it ' • • -J ' 

Pt Pt 

«il r a '^["lL , hi r a s J r [h 1 y , 

Pt Pt 

x Jt r r a s^[fiT >•••]: xg T ~ —a s F[g\T,. . -J , 

Pt Pt 

x /t ~ as-^t/ir 15 ) • • •] ' x 9t ~ "a^^r. • • •] > 

x/i L ~ -2-a s J r [/i^ 1 -',...] , x/t ~ ^-a s J^[/if (1) ,...] , 

xe f ~ — =- a s Fl- • -1 , :re ~ — =- a s T\xe, . . .1 (5.45) 

Pt Pt 

from the 1/p 3 part of <J>(x,pt)- The distributions on the l.h.s. depend on x and pr, and 
on the right-hand side we have convolutions of the form 

F[f] =K q ®f a + K g ®f, (5.46) 

where a is the quark or antiquark flavor of the pT-dependent functions on the l.h.s. of ( 5.44 ) 
or ( 5.45 ). We note that at higher orders in a s one has instead a sum over all quark and 
antiquark flavors, as in (|3.8| ). The contribution from gluon distributions in ( |5.46| ) is absent 
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for chiral-odd distributions. The kernels K q and K g are of course not the same for the 



different functions in ( |5.44 ) and ( 5.45| ), and we use J- in a generic sense. We have explicitly 
calculated the hard-scattering kernels H2 and H3 defined by ( 5.3C| ) and verified that they 
give nonzero contributions for the functions given as arguments of T in (5.44) and ( 5.45| ). 
By ... we have denoted contributions from three-parton correlation functions. We have 
not listed the twist-three distributions parameterizing $3(2;) as arguments of T, because 
they can (up to quark mass suppressed terms) be related to the functions in $9,3 (x) and 
to quark- gluon-quark correlation functions, as we remarked after Eq. ( |5.20 ). An exception 
is e(x), which has no counterpart in $^3(2;). Furthermore, there is no distribution in $2 
or &o,3 that is multiplied by Sl and both T-odd and chiral-odd. At order a s and 1/p 3 the 
high-pr behavior of e L (x,p T ) can hence only be generated from the <3?a,3 term in Q(x,pt)- 
The relations in Q5.44 ) only involve collinear functions of twist two and those in ( |5.45| ) 
only collinear functions of twist three, corresponding to the respective order in the 1/p 
expansion of the correlation function $>(x,pr)- On the other hand, there are pr-dependent 
functions of twist two and twist three in both ( |5.44 ) and ( |5.45| ). In other words, the twist 
of the collinear distributions and the pr-dependent distributions in the high-p^ limit need 
not be the same. 

By power counting, the pr-dependent distributions f^, g- 1 , e T , can receive a con- 
tribution from the 1/p 2 part of $>(x,pr), but explicit calculation at order a s gives a zero 



result. This readily follows from our discussion above ( 5.44 ) because these distributions 
are T-odd. Their high-p-p behavior therefore starts at order a 2 /p 2 and reads 



x 1l ~ ~jo? s T[g^[ , 
Pt 



1 2 
xe T ~ — — a 



Pt 



xg 1 ~ -Y^f/J , 
Pt 

xe? ~ — a 2 F\h-i\ ■ 
Pt 



(5.47) 



There will also be contributions to these functions from <&(x,pr) order a s /p 4 , where the 
necessary T-odd effects can come from lowest-order graphs with three- or four-parton cor- 
relation functions. By power counting in 1/p, these contributions are subleading compared 
with ones in ( |5.47 ), although they appear at lower order in the a s expansion. We omit 
them in our subsequent discussion, but they can easily be restored. 

Explicit calculation also reveals that neither hjj,(x,pr) nor the combination Ii t (x,pt)+ 
hjt(x,pr) receives contributions at order a s and 1/p 2 , although this would be allowed by 
power counting. We shall not investigate the reason of this here, and simply write 
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M 2 . 

r a; 

P T 



1 



P T 



-a 2 

2 s 



r[h x ] . 



(5.48) 



Again there will also be power suppressed contributions at lower order in a s , which go 
like a s /p® for h XT and like a s /p 4 for h T + hjt. We caution that without a full calculation 
of the graphs with multi-parton correlators we cannot exclude that the contributions to 
the distributions given in ( 5.45| ) vanish when all terms are added up. A corresponding 
caveat applies to the a 2 contributions in (|5.47|) and ( 5.48 ). For $>(x,pr) at order a s and 
1/p 2 we give complete and explicit results in section ^. As for the 1/p 3 part of Q(x,pt), 
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the explicit calculation in [6,8] gives f^p{x : p^) 
denotes two- variable convolutions of the form 



(M 2 /p^,)a a J : [G F ,G F ], where T now 



/ 



dxi dx 2 - \r< ( x x 

— — K(xi,x 2 )G F [ — ,— 

Xi X 2 \Xi X 2 



(5.49) 



Given that f^p (x) = G F (x,x), the structure of our result for f^p{x,p\) in ( 5.45| ) 
is hence consistent with the full calculation. 3 

At this point we briefly return to the question of ultraviolet divergences in collinear 
correlation functions, which we mentioned after ( 5.1 2| ) . With the high-p^ behavior given 
m ( ET44D and flOgD one explicitly sees that the integral J d 2 p T $>(x,pr) diverges loga- 
rithmically at high pr, both for the twist- two and the twist-three part of §(x,pt)- The 
corresponding ultraviolet subtractions in the collinear correlator <&(x) result in a logarith- 
mic dependence on the subtraction scale \x for all distributions in ( 5.1 3| ) . This dependence 
is described by DGLAP equations, whose evolution kernels are closely related to the ker- 
nels appearing in the convolutions of ( 5.44| ) and ( [5 .45 ). With ( p. 45 ) one also finds that 
the integral J d 2 p T p^<^ 2 (x,pT) diverges logarithmically. This leads to DGLAP equations 
for the p\ moments of twist-two distributions in the parameterization ( |5.19D of <&q(x), 
which have been investigated in [57]. In contrast, the integral f d 2 p T p^^ 3 (x,pT) diverges 
quadratically in p F according to ( |5.44| ) and ( [5.47 ). In a proper definition for p^-moments 
of twist-three distributions, such as /^^(x) or g^'(x), one would hence have to deal with 
power- like divergences. 

In the dimensional analysis following ( |5.42| ) we have ignored that the hard-scattering 
kernels also depend on the regularization parameter which is Lorentz invariant and has 
mass dimension two. In applications of low-q F factorization one needs £ comparable to the 



large scale, as already mentioned in section 3.1, so that we can restrict our attention to 
( 3> p\. Terms in the hard-scattering kernels going with a positive power oip\jC, are then 
negligible. In contrast, terms with a positive power of C/Pt wou ld lead to a faster p\ falloff 
than derived in this section. They would correspond to power-like rapidity divergences in 
$>(x,pt)- In the explicit calculations at order 1/p 2 in section |8] we will not encounter such 



terms, obtaining only a modification of the power-laws in ( 5.44 ) by logarithms ln(£/pj,). 
A corresponding statement holds for the calculation of f^j>{x,p?n) in [6,8]. 

The high-fey behavior of the fragmentation correlator A(z, k F ) can be obtained in 
full analogy to the case of &{x,pt)- One can readily obtain results by crossing the hard- 
scattering graphs calculated for the distribution functions, replacing x — > 1/z andp F — > hr- 
This gives 



(1) 



D 1 - 1 

~ a s F[D x ] 

An r\jrp 
6 t\jrp 



The relation between f^ 1 ^ and G F can be obtained by combining (5.19) in the present work with 



Eq. (2) in [43] and Eqs. (29), (40) in [13]. Corresponding relations using different parameterizations have 
been given in [56] and [5]. 
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(5.50) 



Compared with their analogs ( 5.46| ), the convolutions 



m 



i 



D a ®K a + D 9 ® K Q 



(5.51) 



have an additional factor 1/z 2 , which reflects the corresponding factor in ( |5.21 ) . 
5.4 Results for structure functions 

Let us begin this section by recalling the expressions for SIDIS structure functions at low 
qr in terms of transverse-momentum-dependent distribution and fragmentation functions. 
Extending earlier work in [11, 12], the study [14] has given a complete set of results at 
leading and first subleading order in 1/Q, i.e., at twist-two and twist-three accuracy. A 
detailed investigation of color gauge invariance and the appropriate choice of gauge links 
has been given in [13]. The calculations just quoted take into account tree- level graphs, 
where gluons are restricted to be collinear to the target or to the observed hadron h and 
only appear when they are attached to the distribution or fragmentation correlators (see 
Fig. 2 in [14]). 

For a compact presentation of the results, we introduce the unit vector h = —Qt/WtI 
and the transverse-momentum convolution 

C[wfD] =J2 xe lJ d 2 p T d 2 k T 5^{p T -k T + q T )w(p T ,k T )f a (x, P 2 T )D a (z,k 2 T ), (5.52) 



where w(p T ,k T ) is an arbitrary function and the sum runs over quarks and antiquarks. 
The results for the structure functions appearing in ( |2.3| ) then read [14] 
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(5.53) 
(5.54) 
(5.55) 
(5.56) 
(5.57) 
, (5.58) 
(5.59) 
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F LL =C[g 1L D 1 ], 
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(5.60) 
, (5.61) 
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(5.62) 
(5.63) 
(5.64) 
(5.65) 



(5.66) 



(5.67) 
(5.68) 



(5.69) 



(5.70) 



In the entries for F uv - L and F^pf^ ^ we have indicated that these structure functions 
come out to be zero when the calculation includes only terms up to order l/Q. The 
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fragmentation functions with a tilde are given by 
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(5.72) 
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= ~~w h D ^ 


(5.73) 
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H k 2 


(5.74) 



Using ( 5.50| ) and neglecting the small contributions proportional to the quark mass m, we 
readily see that the behavior for kx S> M is the same for corresponding functions with and 
without a tilde. 

The tree-level calculations in [11,13,14] do not take into account soft gluon exchange or 
virtual corrections involving hard loops, so that the soft and hard factors we encountered 
in ( |3.3j ) and ( |3.22 ) do not appear in the convolution ( 5.52j ). Detailed investigations of 
factorization for SIDIS with measured qr have recently been given in [26,36] and [27,50], 
extending the seminal work of Collins and Soper [24] . The factorization formulae discussed 
in these papers have the form ( |3.22| ) and are valid at all orders in a s but restricted to the 
leading order in 1/Q. Although a number of subtle issues remain to be fully clarified [27] , we 



will use ( 3.22 ) in the following. Since we aim at deriving expressions at lowest nonvanishing 
order in a s , we can neglect the hard factor \H\ 2 = 1 + 0(a s ). The convolution in ( 5.52; ) 
should then be extended to 

C [wfD] = s £^xe\ I d 2 p T d 2 k T d 2 l T 5 {2) (p T - k T + l T + q T ) 

x w(p T , k T ) f a (x,p 2 T ) D a (z, k 2 T ) U{l 2 T ) . (5.75) 
At high transverse momentum It S> M the soft factor behaves as U(l T ) ~ a s /Z^, with a 



coefficient we shall give in (8.51) below. Our normalization convention is 

d 2 l T U{l' T ) = l + 0(a s ), 



(5.76) 



where it is understood that the integral must be suitably regularized at large It- 

Whether Collins-Soper factorization can be extended to structure functions that are 
of order 1/Q is not known. We note that the study of color gauge invariance in [13] 
was limited to gT-hitegrated observables in this case, and that a problem with twist-three 
factorization has been found in a spectator model calculation [15]. In the following we 
adopt the working hypothesis that the twist-two factorization formula can simply be taken 
over at twist-three accuracy, using the convolution (5.75) also for evaluating the high-g^ 
behavior of the 1/Q suppressed structure functions in ( 5.53 ) to (|5.70|) , We will return to 
this point at the end of section |3,3l 

We now show how to calculate the high-g^ behavior of the convolution ( |5.75 ). At 
order a s , only one of the factors f(x,p T ), D(z, k T ), U(l T ) can be taken at high transverse 
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momentum. Let us first consider the simple case where w(p T , k T ) = 1. In the region where 
Pt is large, we use the 5 function in ( |5.75 ) to perform the px integral and approximate 
p T = k T — l T — q T ~ — q T in f(x,p T ). The remaining integrals over k,T and It can then be 
carried out independently. According to ( |5.21 ) and our discussion after (5.16), the integral 
over kT gives a collinear fragmentation function, up to a s -corrections that can be neglected 
to our accuracy. Likewise, the integral over It gives unity up to a s -corrections according 



to (5.76). Since we are considering the region where Ut and It are small compared with qT 
the integrals over these momenta should be suitably cut off, as is required for ( |5.21 ) and 



(|5.76[) to make sense. Repeating these arguments for the cases where Ut or It are large, 
we obtain 



/ d 2 pT d 2 k T d 2 l T 5^ {pt — k T + l T + q T ) f(x, p T ) D( )U{1 
« /(*> q T ) ^ + f(x) D(z, q T ) + f(x) ^ U{q 2 T ) . 



(5.77) 



For nontrivial functions w(p T ,k T ) the calculation is slightly more involved. Instead of 
approximating e.g. p T = k T — l T — q T ~ — q T , we need to Taylor expand the functions of 
p T around — q T . We take as an example the convolution C[[k T ■ p T ) h^H^ appearing in 
FljU h and consider the region where pr is large. We perform the integral over pt using 
the 5 function and obtain 
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(5.78) 



where both terms in square brackets behave as l/q T . The . . . represent contributions from 
the regions where kT or It is large, which are of the same order in 1/qT- 

As we see in (jU7D , (|5l8|) , and Q ), the leading power behavior of some distribution 
or fragmentation functions comes with a factor o^. At this order, one must also take into 
account regions of integration in ( |5.75| ) where two out of the three momenta pt, kT, It ar e 
large, but it turns out that these do not contribute to the a 2 terms given in the following. 
Using the high-transverse-momentum behavior in (5.44) to (5.48) and ( 5.50| ), we obtain 
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(5.79) 
(5.80) 
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(5.81) 
(5.82) 
(5.83) 
(5.84) 
(5.85) 
(5.86) 
(5.87) 
(5.88) 
(5.89) 
(5.90) 
(5.91) 
(5.92) 
(5.93) 
(5.94) 



Here either the parton distributions or the fragmentation functions are convoluted with 
kernels Ki or Li : 



where the sum runs over quarks and antiquarks for a and over quarks, antiquarks and 
gluons for i. As we will see in section || these kernels contain logarithms of Q/qT- Their 
origin is the dependence of fi(x,p T ) or D\{z, k T ) on ( or (h, which we tacitly omitted 
in ([5.75 ). When resummed to all orders in a s in the way we sketched in section ||, these 
logarithms can lead to a substantial modification of the power laws in (f[7|) to flsM) . 
A numerical study of these effects on azimuthal asymmetries in Drell-Yan production has 
been performed in [58]. 

We note that for the 1/Q suppressed structure functions in ( |5.79 ) to ( 5.94| ), contribu- 
tions from U (l T ) taken at It ~ — Qt ar e power suppressed or have the same power behavior 
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as contributions where either pt ~ — qr or fey ~ qx- For these structure functions, the 
power behavior at high qx hence remains the same if we simply ignore the soft factor and 
work with the tree-level convolution ( |5,52| ) instead of ( |5.75 ). 

6. Comparing results at intermediate q T 

We can now compare the results for the region M <C qr "C Q obtained in the low-gr 
calculation of the previous section with those obtained in the high-qy calculation. As we 
mentioned in section not all structure functions have been calculated in the high-g^ 
picture. For the cases where results are available, we find 

1 

f uu,t ~ — oisF[hDi] , (6.1 
<?t 

F U u,L~-j^a s F[f l D 1 ], (6.2 

Fuu^-^^^ifiD,], (6.4 

nu^-^^TihD,], (6.5 

F LL ~ ~2 a s ^[diDi] , (6.6 

C^^MJ. (6-7 

pffitH-ts) ^^L as ^[G FD ^G F D x ], (6.8 

Fp^-^La^iGrD,], (6.9 

f&UH+ts) ^^£ as ^ hi E F ] : (6.10 
It 

Fffi**-* B) ^-^a s T[G F D ly G F D x l (6.11 
Q z qr 

F ut* S - -^as^lGpD^GpD^hxEp], (6.12 

jfiWH-is) ^^L^jr^^GpD,]. (6.13 

The symbol J- in ( |6.1| ) to ( |6.7] ) has the same meaning as in ( 5.95| ), whereas for the terms 
involving the three-parton correlation functions Gf and Gf we have 

T[GD] =4 [(K i 0G t )(x) J D a (z) + G a (x,x)(^0L i )(z)] , (6.14) 

where the two- variable convolution [Ki <g> G*) (x) is of the form ( |5.49|) . The terms involving 
the three-parton fragmentation function Ep in fl6.10|) and ( 6.12[) are defined in analogy to 

dm. 
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The results in ( |6.lD to (6.4) and (|6,6j ) to ( |6,7D are directly taken from our expressions 
( 4.23[ ) to ( [4.28 ), whereas the result for Ff^^ h in (S.5) has been extracted from the calcula- 
tion in [45]. The form of the Fut structure functions in ( |6,8|) to ( 6.13| ) can be obtained by 
taking the limit qx <C Q of the results of Eguchi et al. [43,44], with the caveats discussed in 
section 0. We note that the results of [43,44] also contain terms involving the product h\Fp 
in Fy^fj! 1 as well as terms involving GpDi or GpD\ in F^ T ^ h+< ^ s \ However, these 
contributions behave like M/(Q 2 qr) for qr <^ Q and are thus power suppressed compared 
with the terms given in ( |6.8[) and ( |6,10| ), 

Let us first discuss the unpolarized structure functions. Comparing the high-gr results 
^1]) to |[D with the \ow-q T results (§7T9|) to ( |5T8l] ) , we find that at intermediate qx the 



power behavior of both F UUT and F^^ h agrees in the two calculations. We shall see in 



section |8.3| that in the case of Fuut this agreement extends to the explicit expression of the 
structure function at order a s . By contrast, the leading power behavior obtained for F™^ 2 ^ h 
in the intermediate region is not the same in the low- and the high-gr calculations. In fact, 



the two results ( |6.4[ ) and ( 5.81 ) describe two different physical mechanisms, since the low-qx 
calculation involves chiral-odd distribution and fragmentation functions, whereas the high- 
qT calculation involves chiral-even ones. Finally, the longitudinal structure function Fuu,L 
only appears at order 1/Q 2 in the low-qx calculation and is hence beyond the accuracy 
of the results given in section [T4|. We remark that it is far from clear whether small-gr 
factorization still holds at twist-four level, given that even the twist-three case is not fully 
understood. 

At this point we wish to discuss the calculation of the unpolarized structure functions at 
low transverse momentum in the parton model [59], where intrinsic transverse momentum 
is included in distribution and fragmentation functions and the kinematics is taken such 
that the quarks in the parton-level subprocess 7*g — > q are on shell. Using Eqs. (4) and 
(32) of [59] and expanding in 1/Q, we obtain Fuu,T = C[fiD\\ as in ( 5.53 ), whereas up to 
relative corrections in 1/Q the other unpolarized structure functions read 
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(6.15) 



with the tree-level convolution C defined in (5.52). The modulations in cos0/i and cos 2<ph 
obtained in this calculation are often referred to as Cahn effect [2,3]. Taking the limit 



qx 3> M of the expressions in (|6.15[) we find the same power behavior as in the high-g^ 
expressions fl6.2|) to (S.4). However, the high-g^ behavior of ( |6.15| ) comes only from the 
high-pj- tail of f\ but not from the high-A^ tail of D\. It hence only involves terms of 
the form (i<Q ® f[) Df, with the same kernels Ki for F^ 



cos (f> h rpcos 2(f> h 
" 7 



UU 



uu 



and F UUL . This 



readily implies that at intermediate qx the parton-model results ( |6.15| ) do not match with 
the explicit results ( [4.24 ) to ( 4.26 ) of the high-g^ calculation. 



We remark that the high-g^ limit of the full twist-three result ( |5.55 ) for F^^ h comes 
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from the chiral-even terms 



2M 



hk T D hp T 
~~M~ z~ + M 



xf ± D l 



(6.16) 



As observed in [14], this coincides with the parton model result (6.15) if one makes the 
Wandzura-Wilczek approximation, i.e., if one sets to zero the functions D = D 1 - — zD\ 
and xf 1 - = xf ± — /i, which are related to quark-gluon-quark correlation functions by the 
equation of motion for the quark field. We will see in section || that for qx 3> M these 
functions are in fact not negligible compared with D\ and fx, so that the approximations 



leading to ( 6.15 ) are not adequate in this limit. In a similar way, one may understand the 



parton model results for F UUL and p^ 2 ^ h as part of the (unknown) complete twist-four 
expression in a low-g-r calculation. They have the correct power behavior to match the 
results ( [4.24 ) and ( |4.26 ) of the high-g^ calculation, but do not reproduce all terms in these 
results. We note that F??? T 2 ^ h has the form (|1.6|) discussed in the introduction. The term 



with coefficient ^2,4 is given by the low-gr result ( |5.81 ), the term with hi,± by the high-gy 
expression (4.26), whereas the parton-model result (6.15) contributes to the subleading 
term Z4 2 in the low-gy power counting scheme. 

Several phenomenological analyses, for instance those in [60-64], have used the parton 
model expressions for the unpolarized structure functions together with the high-g^ results 
( 4.23 ) to ( 4.26|) . We point out that in these papers a Gaussian behavior fi{x,pj) °c 

is assumed for the distribution and fragmentation 



exp[ 



and Di(z, kj,) oc exp[- 



functions appearing in the parton model calculation. Such an approach differs from the 
one taken in the present work, where the power-law behavior of fi(x,p^) and Di(z, kj,) at 
large transverse momentum is retained and explicitly calculated using perturbation theory. 

Turning now to polarized observables, we find that the structure functions F LL and 
pcosiph kave the same power behavior in the high- and low-gy calculations, as do their 



LL 

unpolarized counterparts F UUT and F^^ h . 



As in the unpolarized case, we will see in 
section 8.3 that F LL matches exactly at order a s in the two calculations. Our low-g^ result 
( 5.82[ ) for the T-odd structure function F^^ s has the correct structure to match the limit 
( ]6.5| ) of the calculation at high gy and order a 2 , in [45]. One may expect that our low-gy 
result ( |5.83 ) for F^^ s would also match with a high-gr calculation at the same order, but 
were are not aware of such a calculation in the literature. 

For transverse polarization observables we compare Eqs. (|6l| to (l5l3|) with (pTgTp to 
d-5.91 ) and see that four out of six structure functions have a matching power behavior, 
namely F^ h ~^\ F s J^ hUs) , F S J^ S , and i^ Wh_ * s) . The distribution and fragmenta- 
tion functions appearing in the respective results are compatible as well, given that /{y 1 ^ is 
related with Gf and H^ 1 ^ with Ep. As already mentioned in section [D], the explicit low- 
and high-gy calculations of Fy^% h ^ i n [7, 8] found exact matching at order a s for this 
observable. Looking at the remaining two Fjjt structure functions, we see that F^^-fj 1 ^ 
is beyond the accuracy of the low-transverse-momentum results. This is just as for Fjju,l, 



which is the only other structure function in (2^3) that involves purely longitudinal polar- 
ization of the virtual photon [14]. The structure function F^^"^ h ^ does not match in 
the low- and high-g-r calculations. As in the case of F^ 2 , the low-transverse-momentum 
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result involves chiral-odd functions, whereas the high-transverse-momentum expression in- 
volves chiral-even ones. The low-q T result ( p9|) for F^ 6<ph 0sj can potentially match a 



high-transverse-momentum calculation at twist three and order a 2 , and the high-g^ result 
Q6.ll ) could match with a low-gr calculation at twist four. Both types of calculation are 



beyond the current state of the art. 

To the best of our knowledge, F^ h ~^ s) j p^4>s and F ^(^>h-4>s) have not been 



computed in the high-g^ approach. From the overall factor M in ( |5.92| ) to ( |5.94 ) we can 



only conclude that these low-g^ results can potentially match with those of a high-g^ 
calculation at twist-three accuracy. 

In table Q we collect the results for the leading power behavior of all structure functions 
we have discussed. We notice that for several observables the twist of the low-qx and the 
high-(/T calculation is not the same, which is reminiscent of a similar observation we made 



for the high-pr behavior of distribution functions in section 5.2. 



6.1 Interpolating from low to high qx 

Let us now see how one can practically proceed when the leading terms in the low- and high- 
qT descriptions of an observable do not match in the intermediate region. As an example 
we take the unpolarized structure function F^ 2 ^ h . We denote its Iow-^t approximation 



given in ( 5.56| ) by L^ 2 ^ h , and its high-gr approximation ( |4.26| ) by . Since in the 



intermediate region the two expressions describe distinct contributions to the cross section, 
one may consider to use 

TpCos2(f> h T cos2ij> h TjCos24> h / R ^ ~\ 

as an approximation for this observable. The quality of this approximation can be assessed 
from the power behavior of its terms in the different regions: 

Luu* h ~ Qt/M 4 for q T <M, (6.18) 

~ M 2 jq\ for q T » M , (6.19) 

cos 2 
UU 



H™* 2 ^ ~ 1/Q 2 for all q T , (6.20) 



where the behavior in ( 6.18 ) reflects that L™^j 2 ^ h must vanish like q\ for q T — ► due to 
angular momentum conservation [42]. In the intermediate region M <C qT *C Q both terms 
in ( |6.17| ) are required: together they give an approximation with relative corrections of order 
M 2 /q^ or q T /Q 2 ■ The relative weight of the two terms in this region is L™^j 2 ^ h / H^ 2( ^ h ~ 
M 2 /q^ x Q 2 /q T and thus varies from values above to values below 1. As an aside, let us 
comment on the use of a transverse-momentum-dependence like h^x^p 2 ^) oc exp[— cp 2 ^] 
and H±(z, kj,) oc exp[— Ck 2 -], which is often taken in phenomenological analyses. Whereas 
at small transverse momentum a Gaussian behavior of distribution and fragmentation 
functions is found to give a good description of data in many situations, it misses the 
perturbative tails of these functions. As a result it does not give a good approximation 
of F^ 2 ^ h at intermediate qx- For M < gr < yj MQ the contribution ( |6. 19 ) from the 



perturbative tails is actually dominant, and for a/ MQ < qx <C Q it is only suppressed 
compared with ( 6.20| ) by a factor much larger than M 2 jq\. 
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Table 2: Leading power behavior of SIDIS structure functions in the intermediate region M <C 



Qt <C Q, corresponding to the expansions in (1.2) and (1.4), respectively. Empty fields indicate 
that no calculation is available. The specification of twist 4 for F UTJ L and reflects that 

these observables are zero when calculated at twist-two and twist-three accuracy. 



For large qr ~ Q the ansatz (|6.17| ) can be used as well: the low-q? calculation is 
not valid in this region, but the term L^ 2 ^ h is power suppressed by M 2 /Q 2 compared 
with the leading term H^ 2 ^ h , which itself provides an approximation of F^ 2( ^ h up to 
M 2 /Q 2 corrections. Adding L™^j 2 ^ >h in this region hence does not spoil the accuracy of 
the description. Likewise, the high-g^ calculation is not justified for qr ~ M, but in this 
region the term H^ 2 ^ h is suppressed by M 2 /Q 2 compared with the correct approximation 

does not vanish like q T . 



L C yij 2 ^ h ■ However, one cannot use ( 6.17 ) for qr 



since Hjju 



To repair this, one may instead take 



pcos 2<fr h 



L UU + P 



q T 

M 2 



rrCOS 2<f) h 

n uu 



(6.21) 



with an interpolating function p(r) that satisfies p(r) ~ r for r — > and p{r) — 1 



for 
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r > 1. A simple choice is p(r) = r/(l + r), but obviously there are other possibilities. 
An often considered observable is the azimuthal asymmetry 



4 cos 2<j> h 
A UU 



eF, 



,cos 2<j>f, 



UU 



^UU,T + £ ^UU,L 



(6.22) 



Depending on qj> we can approximate its denominator using the high-g^ expression i?[/{/T+ 
eHjju,l from ( [4.23[) and ( 4.24j ) or the low-qx result Ljjjj^t given in ( 5.53| ) . Since Fjjjj^ is 
suppressed by l/Q 2 for qx <C Q, we do not need the unknown low-qx expression for this 
structure function. Using 



1/M 2 
l/q 2 T 

Huu,t + zHuu^l ~ 1/ Qt 



L>uu,t ~ 
Luu,t ~ 



for q T <M, 
for q T > M . 
for all qx 



together with ( gig ) to fT2C| ), we find that 



/.cos 2</> h 



eL' 



cos 2</> h 
[/[/ 



+ 



cos 2<t> h 
UU 



UU,T 



^UU,T + £ ^UU,L 



(6.23) 
(6.24) 
(6.25) 



(6.26) 



gives a good approximation of the asymmetry in the full qx range. In the intermedi- 
ate region, the denominators of the two terms in ( 6.26j ) coincide up to terms of order 
M 2 jq T or q\IQ 2 and approximate Fjjjj^t + £Fuu,L with that precision. As discussed 
above, both the low-qT and the high-g^ contributions are important in the intermediate 
region M <C qT <C Q (where again one finds that with a Gaussian ansatz for the transverse- 
momentum-dependence of distribution and fragmentation functions, the low-g^ term would 
not be correctly described). For qx ~ Q the low-qx term is power suppressed and may 
hence be kept in ( 6,26j ). For qx^M, the high-g^ term in the asymmetry is suppressed 
by a relative factor of M 2 /Q 2 compared with the 1ow-</t term and does not degrade the 
quality of the approximation ( |6.26 ) in the limit qT — > 0. An additional suppression factor 
as in ( |6.21| ) is therefore not required. Recalling the discussion after Eq. ( [4.33 ), we can 



understand why eH^ 2 ^ h /{H UUT + eH UUL ) has the correct qT — > limit required by 
angular momentum conservation: the propagator factors l/q T that lead to an unphysical 
behavior of the individual structure functions cancel in this ratio. 

Let us finally remark that the discussion in this subsection is at the level of power 
counting arguments. When using ( 6,21| ) or ( |6.26| ) in practice, one can explicitly check 
whether the terms that are out of their region of validity (the L terms for qT ~ Q and the 
H terms for qT < M) are numerically small compared with the leading ones. 



7. Integrating over qT 

7.1 Behavior of integrated and weighted observables 

Up to now we have focused on the gy-dependence of the structure functions F{Q,qx). 
As we mentioned in the introduction, observables that are integrated over qx, with or 
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without a weighting factor ((/t/M) p , can be preferable to observables differential in qx for 
experimental reasons. Without dwelling on such practical issues, we now use our results of 
the previous sections for discussing the theoretical interpretation of integrated observables. 
As a shorthand notation we introduce 



irz 



Q 2 

'I m a x 



dq 2 T (^) P F(Q, qT ) 



(7.1) 



where (/ max is the upper kinematic limit of qx, to be treated as a quantity of order Q in the 
power counting. The prefactor has been chosen for later convenience — note that irz 2 dq T 
corresponds to d 2 Ph±- 

To make the notion of "intermediate transverse momentum" more precise, we introduce 
two scales FM 2 and 7Q 2 such that r > 1, 7 < 1, and TM 2 < jQ 2 . In the intermediate 
region TM 2 < q T < 7Q 2 the results of both the low-qx and the high-g^ calculations 



are then valid, and one can use their respective limiting expressions given in sections 5A 
and ||. It is easy to determine the power-law behavior of the contributions from the regions 
q\ < TM 2 and q\ > 7Q 2 to an integrated observable. For a single term in the general 



low-gr and high-gr expansions ( |1.1|) and (1.3), we obtain 
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M 2 
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T7? / d 1' 
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Or 
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Q 2 

4 max 



T79 / d( l 

1Q 2 



qr_ 
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M 
q T 



MS 



, / IT 



M 
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~Q 



(7.2) 
(7.3) 



from straightforward dimensional analysis. The Q-dependence of the integrals can thus 
be established without knowledge of the functions l n {qT/M) and h n (qT/Q): it is directly 
determined by the twist n in the low-qr case (|7.2|) , and by the twist n and the weighting 
power p in the high-q^ case ( |7.3| ). We observe in particular that for p = 0, i.e. without 
weighting, the twist-two terms in both the low- and high-^T calculations give contributions 
to the integral that stay constant for Q — > 00, whereas higher- twist terms die out in that 
limit. For p > the contribution from the high-g^ region is enhanced: a twist-two term in 
the Iow-^t calculation will then only dominate the integral over all qT if for the observable 
in question a sufficient number of terms with low twist in the high-g^ result are zero. 

As a preparation for the discussion of azimuthal and polarization asymmetries let us 
first take a closer look at the familiar structure functions Fyu^ and Fjjui,. With the 
behavior Fjjjj^t ~ 1/^t in ^ ne intermediate region (obtained in both the low- and high-g^ 
calculations), we obtain 

For the integral in the Iow-c/t domain q 2 -, < TM 2 we have the generic power-law behavior 
given in ( f7.2|) with n = 2 and p = 0. Using in addition that Fjju,t ~ V^t a ^ ^ ne u PP er 
end of the integration region, we have 



TM 2 



dq 2 nF l 



(7.4) 



TM 2 



dq T F l 



UU,T 



In 



r r 



(7.5) 
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with some number To ~ 1. Likewise, we can use that Fjju,t ~ 1/<7t a ^ the l° wer end of 
the integration region in the high-g^ domain q\ > 7Q 2 , and get 



2 

dq^FuuT ~ In — (7-6) 

7 Q2 ' 7 



with some 70 ~ 1. In the complete integral {{Fjjjj,t)) all three regions in (|7.4|) to (\t 
thus contribute at leading power in 1/Q, and the dependence on the artificial separation 
parameters T and 7 cancels as it should. We note that, since we are concerned with power 
behavior in this section, we have not taken into account logarithms of Q/qr in the high- or 
low-gr results for Fuu^t, which would modify the logarithms on the r.h.s. of ( |7.4| ) to (ffTq ) . 

To calculate the integrated structure function one must not double count the contribu- 
tions from the 1ow-</t and high-g^ calculations in the intermediate region. Since the results 
of the two calculations coincide there, one may simply switch from one to the other descrip- 
tion at a suitable point, say at q\ = 7Q 2 . We can now make contact with the standard 
description of gr-integrated SIDIS in the collinear factorization framework, where ((-Fj/[/,t)) 
is expressed in terms of the collinear functions fi(x) and Di(z). Let us in this framework 
take [i 2 = 7Q 2 for the factorization scale and consider the Born graph as well as the real 
and virtual a s -corrections, i.e., the one- loop graphs where a gluon either does or does not 
cross the final-state cut. Loosely speaking, the Born term then corresponds to the sum 
of ( [7.4[ ) and (|7.5| ), and the real corrections to ( |7.6|) . The virtual corrections correspond to 
the hard factor \H\ 2 in the Collins-Soper factorization formula ( |3.3D , which we neglected in 
section 5.3 when extracting results at lowest order in a s . The logarithm of 7Q 2 in ( [7.4| ) cor- 
responds to the scale dependence of the collinear distribution and fragmentation functions, 
whereas the 7-dependence in ( [7.6j ) corresponds to an explicit logarithm ln(Q 2 /n 2 ) in the 
real corrections. At the technical level, however, collinear factorization is typically imple- 
mented by using dimensional regularization instead of a transverse-momentum cutoff. The 
real corrections are then integrated down to qT = 0, whereas collinear distribution and frag- 
mentation functions are defined from integrals f d 2 ~ £ p T fi{x,p\) and J d 2 ~ £ k T D\{x, k^) 
over the full transverse-momentum region. Subtractions defined for instance by the MS 
prescription are then performed, which on one hand ensure that there is no double count- 
ing and on the other hand remove terms corresponding to logarithmic divergences in the 
physical limit e — > 0. Since the incoming and outgoing parton momenta are approximated 
as collinear to the associated hadron momenta, the Born term and the virtual corrections 
appear with a factor 5^ (<7y) in the calculation. 

For the longitudinal structure function Fxju,L the situation is quite different. Using 
the same procedure as for Fyu^ we obtain 

dq T F UUjL ~ 7 1 tw ~ 7, (7.7) 

d lT F uu,L ~ 7o - 7 ~ 1 (7.8) 

7Q 2 



from the result ( |6.2[ ) of the high-g^ calculation. For the second step in (7/7) we have as- 
sumed that (rM 2 )/(7<5 2 ) is sufficiently small compared to 1 — otherwise the corresponding 
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integral becomes small simply because its integration region shrinks to zero. After the 
second step in ( |7,8| ), the dependence on 7 no longer explicitly cancels in the sum of the two 
integrals, but this leads to no inconsistency because (fT7|) is negligible compared with ( fflj| ) . 



One easily sees that the parton-model approximation Q6.15 ), whose power behavior agrees 
with the high-qr result in the intermediate region, gives a result suppressed as M 2 /Q 2 
when integrated over the low-qx domain q T < TM 2 . We thus find that the integrated 
structure function {{Fjju,l}) is dominated by large qx ~ Q and can be calculated from the 
high-gr result alone. Moreover, one can integrate this result down to qx = 0, since the 
contribution from q T < TM 2 is power suppressed by M 2 /Q 2 and thus of the same order as 
the accuracy of the result in the high-gx region. Put differently, one can use the high-gx 
result extrapolated to q T < TM 2 instead of the (unknown) \ow-qx result when evaluating 
the integrated longitudinal structure function. This is just what is done in the standard 
calculation using collinear factorization, where the first nonvanishing contribution to this 
observable starts at order a s - The integration over all qx of the high-gx expression for 
Fuu,L is convergent and simply removes the 5 function in the analog of ( f4.4Q . No subtrac- 
tion is necessary, and correspondingly no dependence on the factorization scale fi arises at 
order a s . 

Let us now turn to the structure functions that describe the ^-dependence of the un- 
polarized cross section. In table || we see that the integrated structure function {{F^^ h )) 
is dominated by large transverse momenta qx, whereas the region where the low-qx cal- 
culation is valid contributes only as a power correction of order M/Q. The condition 
(TM 2 )/(jQ 2 ) < 1 implies VTM/Q < ^7 < 1, so that the factor Vf cannot compen- 
sate the suppression by M/Q. An important consequence is that ((F™^ h )) is not a good 
observable to study the transverse-momentum-dependent distribution and fragmentation 
functions appearing in the low-qx result (5.55). An appropriate observable for this purpose 
is the structure function differential in qx- If integration over qx is required by statistics, 
one should impose a suitable upper cutoff on the integral. According to table |3], the depen- 
dence of the integral on this cutoff is not negligible and must hence explicitly be kept in the 
theoretical calculation. Note that in order not to introduce an artificial (^-dependence, the 
cutoff should be imposed on q T , or equivalently on P%j_, but not on a transverse momentum 
w.r.t. the lepton beam axis. 

Integrated observables which are weighted with a suitable power of qx/M have the 
desirable property that the transverse-momentum convolutions ( |5.52| ) in the low-qx results 
factorize into separate integrals over either distribution or fragmentation functions [12,65]. 
With h = —qx/qr one readily finds from ( 5.55| ) that f dq T (q T / M) F^^ h formally fac- 
torizes into terms involving the p T -moments f^ l \x) and hi^'(x) and corresponding k T - 
moments of fragmentation functions. However, this deconvolution only takes place if one 
integrates over all qx up to infinity. This is clearly inadequate because (q T /M)F^°^ h be- 
comes constant for qx ^» M. A reflection of this is the fact that the p^-moment f^ 1 ^ (x) 
involves a quadratic divergence at large px, as we already noted in section 5^2. Moreover, 
we see in table ^ that the contribution from the low-qx region to (( (q T /M) F^^ h )) is power 
suppressed by M 2 /Q 2 compared with the contribution from qx ~ Q, so that this observ- 
able is even less well suited to study small qx than the unweighted structure function. 
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Table 3: Behavior of selected observables integrated over different regions of q 2 ^ . It is assumed that 
r > 1, 7 < 1 and that {TM 2 )/(^jQ 2 ) is sufficiently small compared to 1. In cases where the \ow-qr 
and high-qT calculations do not match in the intermediate region, their respective contributions are 
given in separate rows. The low-qx entry for Fjju,l corresponds to the parton-model approximation 
in (|l5|). 



Conversely, the weighted structure function is a good observable for studying large qT- The 
high-^T expression for F™^ h depends on the same collinear functions f\(x) and D\(z) as 
Fjju,t but involves different hard-scattering kernels, so that F^^ h provides an additional 
observable if one aims, for instance, at separating the fragmentation functions for different 
quark and antiquark flavors and the gluon, or at testing the adequacy of the theoretical 
description. Up to corrections of order M 2 /Q 2 one can evaluate (((q T /M) F^^ h }} from 
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the high-gr result alone, which in addition may be integrated down to qx = 0. One then 
obtains a simple expression, just as in the analogous case of {{F UUL )). The unweighted 
integral (( ^ h )) is less attractive for studying the high-qy result since the contribution 
from the low-t/y region is only suppressed by M/Q. To evaluate that contribution is difficult 
in practice as it contains transverse-momentum-dependent distribution and fragmentation 
functions that are poorly known. If the weighted integral and the differential structure 
function are affected with large experimental uncertainties, one may instead have to con- 
sider the integral of F^^ h with a lower cutoff on qx- This was for instance done in [17] 
and [62,66]. 

As discussed in the previous section, the structure function F^ 2 ^ h receives contribu- 
tions from the low-qx and high-qr calculations which do not match in the intermediate 
region and have distinct dynamical origins, given that they respectively involve chiral- 
odd and chiral-even distribution and fragmentation functions. As we see in table |, both 
mechanisms contribute to the integrated structure function at leading power, with only 
moderate contributions from intermediate qx- For calculating the integrated structure 
function it is appropriate to add the contributions from the two mechanisms. Further- 
more, it is consistent to perform the (^-integral over the entire kinematical region for both 
mechanisms, i.e. without introducing cutoff parameters, given that contributions from re- 
gions where the approaches are not valid (low q? for the high-gy calculation and vice versa) 
are power-suppressed by M 2 /Q 2 . This is similar to the case of the interpolation formula 
Q6.17 ) discussed in the previous subsection, but for the integrated structure function the 



unphysical behavior of the high-g^ result in the limit qx —* does not matter, at least at 
the level of power counting. The weighted structure function (( (q T /M) 2 F^ 2< ^ h )) has been 
proposed for obtaining a Iow-qt result in terms of the moments h^ 1 ' (x) and (z) of the 
Boer-Mulders and the Collins functions, without any convolution of transverse-momentum- 
dependent factors [12]. According to table ||| this observable is, however, dominated by 
qr ~ Q and only sensitive to hi^'(x) and Hi^\z) at the level of M 2 /Q 2 corrections. 
Such contributions are not under control in the integrated observable, because uncalcu- 
lated corrections of the same size appear in the high-gy region as well. Neglecting M 2 /Q 2 
corrections, one can evaluate (( (q T / M) 2 F^ 2 ^ h )) as an integral of the high-q"r expression 
over the full qx domain. In the same way as (((q T /M) F^^ h }), this provides an indepen- 
dent observable sensitive to the twist-two functions fi(x) and D\{z). 

7.2 Polarization dependence 

Among the many observables for polarized SIDIS, the structure functions F^pp 1 ^ arm 
psin(4> h +<j>s) jj aye received particular attention in the recent literature. According to the 



low-qT results ( |5.62| ) and ( 5.64 ), they provide access to the Sivers function f^f in the 



first and to the transversity distribution h 1 and the Collins function in the second 
case [67]. Both structure functions have been found to be of significant size in HERMES 
measurements on a proton target [20]. 

As we see in table ||, the integrated structure function {{F^ T ^ T h ^ s ^)) is dominated 
by the Iow-^t region and can hence be used for extracting information about fypix^p 2 ^). 
The high-gy region is however only suppressed by M/Q, so that it may be of advantage to 
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impose an upper cutoff on the qx integral in such analysis. The weighted structure function 
(((q T /M) FyT% h <?is ' ) )) receives contributions from both high and low qx at leading order 
in M/Q. One can thus compute the weighted integral by switching from one to the other 
formulation at some qx- To achieve a factorization of the transverse-momentum convolution 
in the low-qx expression, one should however integrate it over all qx up to infinity. Since 
(q T / M) Fy T ^ T h behaves as l/q T for q T 3> M, a suitable regularization is required. This 
suggests a procedure akin to the description of ((Fjju,t)) in collinear factorization, which 
we reviewed in the previous subsection. As dimensional regularization preserves rotation 
invariance in the transverse plane, the integral over all qx of the weighted low-qx result 



( |5.62[) turns into the product 

(( (q T /M) F$ T f-^ )) = -2^xe 2 a (x; Q) D^z; Q) (7.9) 

a 

of collinear functions defined in the MS scheme. One can trade fyx^ ( x ) f° r the twist-three 
function Gp(x,x), which appears in the high-qy calculation [7,8,44]. In order to integrate 
the high-g^ result down to qx = 0, one must extend it to 4 — e dimensions and perform the 
necessary MS subtractions. Adding graphs with virtual corrections to the hard-scattering 
subprocess (which give the hard factor \H\ 2 in the Collins-Soper formalism) one will obtain 
a complete NLO result in a s . Such a procedure would be the analog of a standard NLO 
computation for integrated observables within collinear factorization at twist-two level. 
Note that ( |7.9| ) gives a consistent approximation of the weighted structure function at LO 
in a s , in analogy to the well-known tree-level expression ((-FW,t)) = ^2 a xe 2 fi ( x ) ^1 ( z )- 
The factorization scale /i of the functions in ( |7.9| ) has been set to Q in order to avoid large 
logarithms of Q/fJ, appearing in the a s -corrections. To leading order, the logarithmic Q 
dependence of the weighted structure function then follows from the evolution equations 
for D 1 (z) and /-j™ (or). The latter have been investigated in [57]. 

The situation for the structure function p^^ h+( t' s ' ) [ s the same as for F^ T ^ T h ^ since 
the power behavior of these observables coincides in both the low- and high-qr calculations. 
The evaluation of the weighted structure function in collinear factorization gives 

Uq T /M h )F^ + ^}) = 2 J> e 2 h^x-^H^iz-Q) (7.10) 

a 

for the Born term. The A^-moment H^\z) is related to the twist-three fragmentation 
function Ep appearing at order a s . We note that according to the high-g^ results in [43,44] 
the a s -corrections to both ( |7.9| ) and ( 7.10| ) involve each of the twist-three functions Gp, 
Gf, and Ep- 

According to table [|, the integrated structure functions ((F^ T ^ S )) and {{F^ T ^ 2 ^ >h ^ )) 
receive comparable contributions from all regions of qx- Weighting the structure functions 
with (qx/M) 2 one obtains integrals that can be evaluated in the high-gy formalism up 
to corrections of order M 2 /Q 2 , similarly to the case of (( (q T /M) F™^ h )) we discussed 
in the previous subsection. The high-gy expressions computed in [43, 44] imply that 
(((q T /M) 2 F s J T {2(ph -^ s) )) is sensitive to G F and G F , whereas (((q T /M) 2 F^ 3 )) also de- 
pends on Ep. Whether these observables are large enough to be measured in practice is, 
of course, a different question. 
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From (|5~89|) we can infer that the integral of 



sin(34> h -<f> s ) 
UT 



receives a contribution of 



order 1 from low qj>, whereas the high-g^ result (|6.11|) is suppressed by M/Q. According 



to ( 5.65 ) the integrated structure function \\^ UT 



psm(34> h 4> s ) ^ hence be used to extract 



information on hj- T (x,pT) and on the Collins fragmentation function, with the same caveat 
we discussed for ((F^^p 1 < ^ s ' 1 )). To obtain an integral that is dominated by the high-g^ 



result up to M/Q corrections, one must weight F UT 



,sin(3<j) h -<t>s) 



with (qr/MY 



Let us now turn to observables that involve longitudinal polarization. Similarly to 
Fuu^ h , the lepton-helicity-dependent structure function F^^ h for an unpolarized target 
receives a contribution of order M/Q from low qT and of order unity from high qx- It is 
therefore in principle suitable for investigating the high-g^ result of Hagiwara et al. [45] . 
However, the contribution from large qx comes with a factor a 2 s in this case, which may not 
be sufficient for neglecting power-suppressed contributions from low qx in practice. From 
this point of view, it would be advantageous to weight the structure function with qx/M, 
or to integrate over qx starting from a lower cutoff. 

Finally, the structure functions F LL and F?f ^ h have the same power behavior as their 
unpolarized counterparts F UUT and F^^ h , and their discussion is analogous to the one 
in the previous subsection. In particular, the weighted integral (( (q T /M)F'j°^ h )) depends 
on the polarized parton densities g\ and, if measurable with sufficient accuracy, could be 
used in addition to the well-known observable {{F LL )) for disentangling the contributions 
from different quark and antiquark flavors and from the gluon. 



8. Prom low to intermediate q T : explicit calculation 

In this section we compute the high-transverse-momentum tails of the quark distributions 
in ( [5.44 ) and of the analogous fragmentation functions. These are the functions which 
appear at lowest order in the 1/pt expansion of section [|]^ and are hence expressed in 
terms of collinear functions of twist two. While in section ^ we identified observables 
whose power behavior agrees in the low- and high-g-p calculations, we will then be able 
to check for selected structure functions whether agreement is also found for their explicit 
expressions. 



8.1 High-pr tails of distribution functions 

Let us begin with the quark distribution functions. We work in the original scheme of 
Collins and Soper [24], using a spacelike axial gauge with the singularities of the gluon 
propagator regulated by the principal value prescription. The only Feynman diagrams 
to be evaluated are then those depicted in Fig. |5|a and b. For further discussion and a 
comparison with the calculation in Feynman gauge, we refer to appendix [5} 
The contribution of the quark-to-quark term shown in Fig. |B|a reads 

- ' r dp- I dl + 5{(l - pf) 6(1 



,X,PT) 



§0 (2vr) : 



C F 



x d^(l-p;v) 



x\ p 

~ J 7m "a 

xJ p* 



p 



i-=o,i T =o T 



(8.1) 
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(a) (b) 

Figure 5: Diagrams for the calculation of the leading high-py behavior of the quark-quark 
correlator <&(x,pt) in axial gauge A ■ v = 0. 



where it is understood that p + = xP + and l + = p + /x. As explained in section |5.3| , the 
restriction to leading order in 1/pr allows us to set l~ and to zero when calculating the 
hard-scattering subprocess, and to retain only the twist-two part Q\(x/x) of the collinear 
quark-quark correlator at the bottom of the graph. The gluon polarization sum in A ■ v = 
gauge is given by 

d^(q; v) = -<T + Q q ~ -r 1 ^ v 2 , (8.2) 

q-v \q-vy 

where the singularities at q ■ v = are to be regulated by the principal value prescription. 
Using the (^-function to perform the p~ integration, 

we obtain 

*'^\&r^Y CF ¥ T ft tl - *) J(^»!(|)t,J», (8.4) 

where we have introduced the notation 

2 

p = p + n + - ——^n.+pT, (8.5) 
- V + 

l = i—n + (8.6) 

x 

for the approximated momenta in the hard-scattering kernel. We note that the virtuality 
of the upper quark legs 

f = "T^- (8.7) 
1 — x 

is always spacelike. The gauge fixing vector can be written as 

2(P+) 2 v- , s 

v = v ~ n _ - _i — 1 n+ (8.8) 

with 

( __QP£. = _ *£pr t (8 . 9) 
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where in the second step we have neglected M 2 compared with £. We therefore have 



2{p + ) 2 

(I - p) -v = (I - p) + v 7T— (l—p)~v~ 

x A Q 

where we have introduced the parameter 



1 — x — r] 



1 — x 



p v 



f] 



Pt 



Pt v 



x 2 ( 2{p + ) 2 v~ 
We now decompose the gluon polarization sum as 



3.10) 



3.11) 



dT{l -p;v) =^2d^{l -p;v) , 



i=i 



with 



dty(l-p;v) = j] 
d^(l-P;v) = - 



1 - X 



x 2 g- 



d^(l-p; ra_ 



l-xl-x-rj^-/ 

2xp+ a-pYnX + (I-pYnH 
77— ■ 



Pt 



4x 2 (i-pY(T-p) 



Pt [ 1 — x — rj j 



.12) 



3.13) 
3.14) 
3.15) 
3.16) 



Notice that the first term ( 8.13j ) is proportional to the polarization sum dY v {], — p;n_) 
one would use when calculating in light-cone gauge A ■ n_ =0. We will see shortly that 
the prefactor in ( |S.13| ) regulates the divergence at x = 1 which would arise in that gauge. 
From the parameterization ( |5.13D we readily see that the twist-two part of the quark-quark 
correlator satisfies = ^2 ^+ = 0, so that terms with or n v + in dP" v vanish when 

inserted into (^J). With I being proportional to n + , we hence need only the first two terms 
and the p il p v part of the last term in the decomposition ( |8.12| ). This gives 



<5> q (x,p T ) 



§0 " (2^ ^ / 



1 

P T 
1 — x 



1 dx 
x 



' (l-x) 2 ^(/"-p;n_) #7^I /; '' 



(1 — x) 2 — T)X 2 p?^ 

(1 — x) A — r\x A p T \x 



7m ; 



+ 



77(1 — x) 



[(1 — x) 2 — 7]X 2 ] 2 



.17) 
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To proceed we must determine the behavior of the different terms in the limit x — * 1, where 
p- (i _ becomes singular. Using the form ( |8.2[ ) with v replaced by n_, we obtain 

2 

= -(1-x) 2 p 7 ^h^ + ^f{P^ q 2 + ^ 2 ^p). (8.18) 

Since the minus-component of p drops out in f 4- and 4- Pi the expression in ( |3.18| ) is 
finite for x — ► 1. For the second term in ( 8.17| ) we have 



1 - x ) 1 7 M *I = "(I {ft 4+ - g\S Ll5 4+) f (8.19) 



after plugging in the parameterization of §\ from ( 5,13; ), This contains a piece with two 



factors of p~ , which is proportional to the Dirac matrices 4- or 75 4— ■ According to the 
decomposition fl5,4| ) it therefore does not contribute to the twist-two or twist-three parts 
of the correlator &(x,pt), on which we concentrate here. In the twist-four part of ( |S.1| ) 
this piece leads to a singularity at x = 1, or in other words at p~ — » — oo, showing that 
at twist-four level the A ■ v = gauge is insufficient to render the integral over p~ in 
$(x,pr) = f dp" <&(p) well defined. 

In the following we take the limit p 2 -, -C C corresponding to rj -C 1. The motivation 
for this is that in physical processes we need the correlator <£(x,pr;C) f° r Vr ^ Q 2 an( ^ 
x 2 Q ~ Q 2 , as discussed in section |3.1| . We note that in a frame where xP + ~ Q this 
corresponds to v + ~ v ~. According to ( |8.11 ) the parameter y/rj is then proportional to the 



small angle between the quark momentum p and the hadron momentum P, with a factor 
of proportionality of order 1. Notice that at this point we introduce a hierarchy in size 



between px and P + , which were not distinguished in the power counting of section ^3 
This is similar to what we have done with the high-gy calculation of structure functions in 
section |] : we started with the result (|4.4j) , which is derived without making a distinction 
between the size of qx and Q, and in a second step we took its limit for <C Q. 

For rj <C 1 the first term in the square brackets of ( |3.17| ) can be rewritten by using 
that for any function G(x) which is regular at x = 1 

lim PV f 1 dx 1 ~ X G(x) = C dx G ( X ). + - G(l) In , (8.20) 



77-0 J x (l-x) 2 -rjx 2 J x (l-x) + 2 \rj 

where the plus-distribution is defined as in (P~22j ). From 



lim PV [ dx ^ G(x) = (8.21) 

v-+o J x (1 - x) 2 - r\x 2 

we see that the second term in ( ^.17| ) does not contribute in the small- r\ limit when restricted 
to the twist-two and twist-three parts of $(x,pr)- In contrast, the third term in ( 8.17| ) 
does contribute, since 

lim PV f 1 dx V ^r X \ 9 , 9 G(x) = -- G(l) . (8.22) 



77—0 J x ' [(1 - x) 2 - f]X 2 ] 2 2 
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Therefore, our final result reads 



$ q (x,p T ) 



- Cf— f 1 — 

i) 2vr 2 F p\ h x 

x' 2 



+ -5(1 

2 v 



x) In — 

rj 



1 



x 



2p; 



d""(i-p;n_) ^ lu m-U^-5{l-x) & 2 {- 



.23) 



to leading order in 1/pt, where it is understood that we have restricted ourselves to the 
twist-two and twist-three parts of the correlator on the l.h.s. We note at this point that 
if we work with a timelike axial gauge, i.e. with negative £ and rj in ( |8.8|) to ( 8.11 ), we 



obtain the same result as in (8.23) with ln(— r? -1 ) instead of ln(?7 _1 ). The polarization sum 
d^ w (l — p;v) is then nonsingular in the whole region x < x < 1, and the principal value 
prescription in ( |8.20| ) to (8.22) is not required. A timelike vector v was indeed used for 
the construction of factorization by Ji et al. [26], whereas arguments in favor of taking v 
spacelike were given by Collins and Metz in [50]. 

The gluon-to-quark contribution to the correlation function comes from the diagram 
in Fig. ||b. Its calculation is simpler than the previous one, due to the absence of a gluon 
polarization sum in axial gauge. Correspondingly, the result is independent of rj. The 
counterpart of the expression in (|8.4|) now reads 



<5> q (x,p T ) 



a« „ 1 f 1 dx 

x 



(§b) (27T)2-"p4 



(l-x)^g)| > (/-^) 7 ^, (8.24) 



where the twist-two part of the collinear gluon correlation function is given by 

1 



/IV I 



2xP+ 



.25) 



see e.g. [68]. Inserting ( p.5|) and ( p.6[) and using some Dirac algebra, one finds that the 
integrand of ( |S.24D is finite at x = 1. 

From ( pOaD and (|3|) we can easily project out the contributions to the individual 
terms in the decomposition ( |5.4j ) of ^ q (x,pt)- For the high-py behavior of the unpolarized 
distributions we obtain 



a.s _J_ 
Pt 



2vr 2 - 2 



xf^(x,p 2 T ) = ^ 1 



2 

Lir,- 1 



2tt 2 2p 2 T 



whereas for the polarized distributions we find 



fl{x) - C F fl{x) + (P qq ® fl + P qg ® /f)(a 
^/f(x) + (P 9 >/f + P g >/f)(x) 



3.26) 
1.27) 



a, 1 



xgt q {x,p T ) = 



a, 1 



2 



2vr 2 2p 2 T 



g\(x) - C F g\(x) + (AP qq ®g\ + AP qg ® g{) (z) 



.28) 
.29) 
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(a) 



(b) 



Figure 6: Diagrams for the calculation of the leading high-fcr behavior of the quark-quark 
fragmentation correlator A(z, kx) in axial gauge A ■ v = 0. 



and 



h\{x,p\) 



° Pt 



2vr 2 - 2 



xh^ q {x,p z T ) = — ^ 



1 



2tt 2 2p 2 



2 



ft? (a) - C F ft?(x) + (5P W ® ft?) (a;) 
U?(a;) + (<5P OT ®ft?)(x) 



X ll^-i (^X ^ Prp^j X hrp^ [X ) Pj-t^) , 



5.30) 
5.31) 
5.32) 



where L(r/ x ) is defined as in ( 3.13 ). Here 



<5P w (s) = C F 



L(r^ + 2 5(1 " x) 



5.33) 



is the leading-order DGLAP splitting function for the transversity distribution [69], and 



the remaining splitting functions are given in ( |4.29| ) to ( 4.33; ). The chiral-odd quark distri- 
butions in ( |3.30 ) to ( |8.32| ) receive no contribution from ( 8.24]) because chirality is conserved 
for the quark line in the graph of Fig. [5|b. 

The diagrams for the high-p^ behavior of the antiquark correlation function ^(x,px) 
are obtained from those in Fig. [| by reversing the direction of the fermion lines. The 
results have the form of (8.26) to ( p. 32 ), with identical splitting functions and with all 
quark distributions replaced by antiquark distributions. 



8.2 High-£>r tails of fragmentation functions 

The calculation of the high-transverse-momentum tails of quark fragmentation functions 
proceeds in close analogy to the case of distribution functions. We nevertheless present the 
essential steps in this subsection, so as to show that no problems occur when going from a 
spacelike to a timelike situation. 

The Feynman diagrams to be evaluated in A ■ v = gauge are drawn in Fig. ^. We first 
consider the quark-to-quark contribution of Fig. ^a. The analog of the starting expression 
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8.1| ) now reads 

A q (z,k T ) 



(BO (2tt) 3 



CF~ z j dk + j dr 5((k-lf) 



n 



x^-^)|7^Al(|) 7 ^ 



.34) 



/+=o, i T =0 T 

where k~ = P^ / 'z and l~ = zk~ . The factors 1/z and z/z in ( 8.34j ) arise from the 
definitions ( |5.9| ) and ( |5.21 ) of the fragmentation correlators. We perform the fc + -integration 
using 

1 



5{{k-lf 



i+=o t i T =o T 2k~(\ — z 



5 k 



2k-{l-z) 



and obtain 

A q (z,k T ) 
with 



C F - ' 



§0 (2vr)2 z 2 fc 4 



2k~(l - z) 
I = zkTri- . 



n + + k n_ + /c-r , 



The virtuality of the fragmenting quark 



k 



1 - Z 



5.35) 
5.36) 

5.37) 
5.38) 

5.39) 



is always timelike, in contrast to its counterpart p 2 in the distribution correlator. For the 
calculation of the fragmentation correlator, it is useful to write the gauge vector as 



v = v ' n_i 



Ch 



• n_ 



with 



(2P h -v)\ 2(P h -)*v" 



In analogy to ( [slop and (18.11|) we can then write 

(k-i)-v=(k- i)-v + - (k - i) + v^ 



z- 2 Ch 



1-z-Vh 



l-z 



k v + , 



where we have introduced 



Vh 



k\ v~ 



5.40) 
5.41) 

5.42) 
5.43) 



z - 2 Ch 2(k-y v + ■ 

Taking the limit h\ <C Ch and following similar steps as in the previous subsection, we 
obtain 



A q (z,k T ) 
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.44) 
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to leading order in where as in the case of distribution functions the result is restricted 
to the twist-two and twist-three parts of the correlation function on the l.h.s. For the 
quark-to-gluon contribution, the diagram of Fig. |b gives 



Ai(z,k T ) 



O'.s 



(§b) (2vr) : 



C F 



dz, 



1 - z) A 9 ^ v 



fatf- J) -Tut, (8.45) 



where the twist-two part of the gluon fragmentation correlator is parameterized by just 
one function, 

A^(z) = -^g^Df(z), (8.46) 



2p: 



because we consider an unpolarized hadron. With the parameterization ( |5.9| ) of A q (z, kx) 
we obtain the high-A>r behavior 
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2ir 2 z 2 k\ 
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S.47) 



(8.48) 



from (|8.44| ) and (|8.45| ). For Z)- 19 this implies 



5^(^,4: 



2vr 2 2z 2 4 



£>?(*) - 2C FJ D?(z) + (D? ® P' qq + Df ® /^) (z) 



.49) 



according to its definition ( 5.71| ). Analogous results with the same kernels are obtained for 
the antiquark fragmentation functions Df, D^, and D^. 



8.3 Results for structure functions and their consequences 

We are now ready to compute the behavior of the structure functions Fuut, -^jju h '> 
and F L °£ i ^ h at intermediate transverse momentum. For Fuut we start from the low-qx 
result ( |5.53| ), with the convolution defined in ( |5.75| ). Using the expansion ( |5.77|) we have 



Fuu,t = ^2xe 2 a 



ft(x, &) ^ + ft (x) Dt(z, q 2 ) + ft (x) £M U(4) 



z- 



.50) 



for M < < Q. The high-transverse-momentum limits of ft[x,q^) and Df{z,q^) are 
respectively given in ( 8.26| ) and ( p. 47 ). For the corresponding limit of the soft factor one 
obtains 

2 (x s Cf 1 



U(qj 



tt 2 q\ 



.51) 
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from [24], as we show in appendix |C| Given that 2P + P h = zQ 2 /x up to mass corrections, 
the relations Q, (|Tl]) and flOID , dQa|) imply 

^Q 2 



3.52) 



where in the second equation we have set fc r and p T equal to q T , as appropriate for 
evaluating ( 8.5CQ . Putting the above results together, we obtain 




Fuu > T = ~^7^1^ xe " 



z It a 



/f(x) - C F ft{x) + {P qq ® /? + P gs 8) /f)(a 



L (Vh ) n « 



£>?(*) - CpDf(z) + (£)? ® Pqq + D[ ® P gq ) (* 



+ 2C F / 1 a (x)^(z) 

1 Q 



6., 



Q 5 



,/;;' c. c ) l ( ^ ) + (i>? ^ + ® p gq ) (z) 



+ (P qq ®fi+P q9 ®f!)(x) D<t(z) 



3.53) 



which is identical with the result ( [4.23[) of the high-qx calculation. The same agreement 
has been found by Ji et al. [7], who used the low-q^ factorization scheme specified in [26] 
instead of the original Collins-Soper scheme [24]. Note that the terms with Cpff (x) Df(z) 
cancel among the different contributions in (|8.53| ). By virtue of ( |3.52 ) the dependence on 
the gauge parameters r\ and rj^ also cancels, as it should. We remark that we obtain the 
same final result if we take a timelike vector v instead of a spacelike one. Both r] and % 
are then negative, and L(t]~ 1 ) and L{rj^ 1 ) are replaced by L{—rj~ l ) and L{— r]^ 1 ), so that 
they still add up to 2L(Q 2 /q T ). For f ll 

we obtain a result analogous to fDj5), with the 
parton distributions f\ replaced by g± and the convolutions P ® /i by AP ® g\ . This result 
exactly matches the expression ( [1.27 ) obtained in the high-gy calculation. 



We now turn to the structure function Fi 



cos <f> h 



UU 



According to ( 5.45 ) and (|5.50[) the 



terms with and H± in the low-g^ expression (5.55) are power suppressed compared to 
the terms with / and D when qx ^> M. For intermediate qx we therefore have 



UU 



' Q V 



xf La (x,q T ) 



Df(z) 



D ±a (z,q*)1 



3.54) 



at leading power and leading order in a s . In this case there is no leading contribution 
from the soft factor taken at large transverse momentum. Proceeding as we did in ( 5,7§| ), 



one finds that the leading term in the expansion of U(l T ) around It = —Qt gives zero in 



the convolution ( p. 55 ) because it does not depend on a direction in the transverse plane, 
whereas the next terms in the expansion only give contributions that are power suppressed 



compared to those in ( 8.54 ). We therefore obtain the same result ( 8.54 ) if we omit the 
soft factor in the transverse-momentum convolution ( |5.75D . Using the high-transverse- 
momentum limits ( 8.27| ) and ( |8.49| ) of f ±a {x,q T ) and D ±a (z,q T ), we get 
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.55) 



which is noi identical to the high-g-p result ( |4.25| ) because of the extra term 2Cj?/f (x) Df(z) 
in the brackets. The same situation is found for Fff ^ h , with fi replaced by g\ and P' ® /i 
by AP' ® gi 

This disagreement has important consequences. Since the leading terms in the high-g^ 
and the low-gr calculation of F™^ h have the same power behavior for M <C qx <S Q, 
their explicit expressions in that region must agree if both of them are calculated correctly. 
This is clear, since both calculations give the same term of a double expansion in M/qx 
and qx/Qi as given in (]1.2[) and (|l.4|). We have no reason to doubt the validity of the high- 



ly result ( 4.25 ), which comes from a twist-two calculation in collinear factorization. The 
same holds for the high-transverse-momentum behavior of the functions f ±a (x,q T ) and 
D ±a (z,q T ) in (|8.27|) and ( 8.49 ). In contrast, the Iow-qt expression we used for F™^ h is a 



twist-three result, for which no proof of factorization is available. To obtain the expression 
in ( |8.54j ) we have assumed that the tree-level result (|5.55| ) can be generalized by taking 
over the convolution ( 5.75| ) established for the twist-two sector. The comparison of ( 8.55] ) 
with (4.25) implies that this assumption is incorrect. 

Based on our finding, one may speculate how a correct twist-three factorization formula 
will look like if factorization can be established at that level. Simple modification of the 
soft factor U(l T ) can obviously not yield agreement with the high-g-r result since this factor 
does not appear in the limiting expression ( 8.54 ) for the reasons we explained above. The 



situation would be different if the soft factor were dependent on the direction of It, which 
would require it to have a nontrivial structure in either Lorentz or Dirac space (through 
factors l T or f T ). Such a dependence would go beyond the eikonal approximation for the 
coupling of soft gluons to fast partons, which may be necessary at subleading order in 
1/Q. We shall not pursue such speculations here. Clearly, the requirement to match the 



high-<7T result ( 4.25 ) for F^^ h at intermediate qx can be used as a consistency check for 



any framework that extends Collins-Soper factorization to the twist-three sector. 

It is instructive to note that the low- and high-q-r results disagree by a term pro- 
portional to fi(x)Df(z), where neither the distribution nor the fragmentation function 
appears in a convolution over longitudinal momentum fractions. In the calculations of the 
previous subsections, such terms arise from configurations where a gluon has zero plus- or 
minus-momentum. The correct treatment of this phase space region is nontrivial already 
in proofs of factorization at the twist-two level [26,27], so that it is not too surprising that 
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this is where problems occur in the naive extension to twist three which we have explored. 
At this point we return to the issue of transverse-momentum-resummation for Fr? s 



uu 



which we have briefly discussed in section |3.3| . We can now understand why the splitting 
functions P' qq , P' qg: and P' in the high-gy result ( 4.25| ) are different from the usual DGLAP 
kernels. Up to 5-function terms they describe the high-transverse-momentum behavior of 
/-*- and D^~, rather than the one of the more familiar functions f\ and D\. A corresponding 
remark applies to the cos 4> asymmetry in Drell-Yan production investigated in [9] . If a low- 
qx factorization formula for these observables can be established, it should also allow one 
to adapt the original CSS procedure [1] for the resummation of large logarithms ln(Q 2 /q^) 
at next-to-leading logarithmic accuracy and beyond. From this point of view, resummation 
for F^ 2 ^ h and its analogs in Drell-Yan production or e + e~ annihilation appears rather 
daunting since it would require a formulation of low-qx factorization at twist-four level, 



extending the simple parton-model result in ( 6.15| ) and putting it on a rigorous footing. 



9. Summary 

The description of semi-inclusive deep inelastic scattering with measured transverse mo- 
mentum qT involves two theoretical frameworks: at low qT one has a factorized representa- 
tion in terms of transverse-momentum-dependent distribution and fragmentation functions, 
whereas at high q? standard collinear factorization can be used. We have systematically 
analyzed the relation between the two descriptions at intermediate transverse momentum 
M <C qr <C Q, where both are applicable. Depending on the specific observable, the 
leading terms in the two descriptions may or may not coincide. 

Using dimensional analysis and Lorentz invariance, we have derived the general behav- 
ior at high px for all transverse momentum-dependent parton distributions of twist two or 
three. The results, listed in Eqs. ( |5.44D to ( |5.48D , involve the convolution of collinear parton 
distributions with hard-scattering kernels, which in the simplest cases are closely related 
to the well-known DGLAP splitting functions. We have computed these kernels at leading 
order in a s for those cases where the collinear distributions are of leading twist, obtain- 
ing the expressions (|8.26| ) to ( 8. 32; ). With these results and their analogs for transverse- 



momentum-dependent fragmentation functions we could establish in Eqs. (|5\79p to f[9|) 



the power behavior for M <C qT <C Q of all SIDIS structure functions that appear in the 
low-qT description at twist-two or twist-three accuracy, allowing for arbitrary polarization 
of target and beam. 

In the high-^T description at order a s one finds a considerable simplification when 
taking the limit qx <C Q : the expressions of the structure functions then involve a convolu- 
tion of either the distribution or the fragmentation functions with hard-scattering kernels, 
whereas the other function is evaluated at the momentum fraction x ox z fixed by the kine- 
matics of the final state. For observables where the high-g^ and low-qx calculations match, 
these kernels can be identified with the ones describing the high-transverse-momentum be- 
havior of the functions appearing in the low-qx description. In such a situation one can 
use the procedure of Collins, Soper, and Sterman to resum large logarithms of Q 2 /q T to 
all orders in perturbation theory. A prerequisite for this is that the power behavior of 
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the observable in the low- and high-g^ calculations must match. We have compared the 
corresponding powers for a wide range of observables, using our low-qx results ( 5.7S| ) to 
( |5.94 ) and their counterparts ( |6.1| ) to ( |6. 13| ) for those structure functions that have been 
evaluated in the high-g^ formulation. This comparison, compiled in table |2|, is one of the 
main results of our work. 

When the two formulations give the same power law at intermediate qx for a given 
observable, their explicit results must agree exactly because they describe the same term of 
a double expansion in M/qx and qx/Q- This constitutes a nontrivial consistency check for 
both the low- and high-g^ calculations. Confirming earlier results in the literature, we have 
verified that there is such an agreement for the unpolarized structure function Fuut, as 
well as for its analog Fxl for longitudinal beam and target polarization. By contrast, the 
structure function Fjju,l for longitudinal photon polarization only appears at twist four in 
the \ow-qx framework, where a complete result is not available. A simple calculation in the 
parton model gives a power behavior which in the intermediate region matches the one of 
the well-established high-g^ result but fails to reproduce its exact form. 

A more involved picture arises for azimuthal asymmetries, even in unpolarized scat- 
tering. At low qT the structure function F^ 2 ^ h is expressed in terms of the Boer-Mulders 
function h± and the Collins fragmentation function H^, both of which are chiral-odd, 
whereas the high-g^ expression involves the usual unpolarized distribution and fragmen- 
tation functions /i and D%, which are chiral-even. The two results thus describe different 
physical mechanisms, which is consistent with our finding that at intermediate qx they 
have a different power behavior. In this region, the two results may hence be added. In 
practice, some arbitrariness is involved in deciding what "intermediate" qT values are. We 
have shown that the sum of the high-gr and the \ow-qx expressions gives a valid approxi- 
mation for F^ 2 ^ h also at large qx, where the low-qx result cannot be trusted but is power 
suppressed compared with the high-qr expression. The latter, however, fails to vanish in 
the limit qx — ► 0, as required by angular momentum conservation, and should hence not be 
used at low qx- A more favorable observable in this respect is the cos2(/>/i asymmetry, i.e., 
the ratio of F^ 2 ^ h and the <j>h independent part Fjjjj,t + £Fjju,l of the cross section. In 
this case, the sum O) of the expressions calculated for low and high qx gives a consistent 
approximation for all transverse momenta, up to corrections of order M 2 /q T and q T /Q 2 - 
The result of a parton-model calculation at low qx, often referred to as Cahn effect, has the 
same property for F^ 2r ^ h as it has for F uv L : its power behavior agrees with the high-g^ 
result in the intermediate region, but its explicit expression does not. The parton-model 
result may hence only be regarded as a partial estimate for the full but unknown twist-four 
correction to F^ 2 ^ h at low qx- 

The description of the structure function F^^ h is more problematic: at high qx it can 
be evaluated in collinear factorization at twist-two level, but at low qx it requires a twist- 
three calculation, for which transverse-momentum-dependent factorization at all orders in 
a s has not been established. As a working hypothesis we have taken the well-established 
result of a tree-level calculation at low qx and assumed that the soft factor which explicitly 
appears in the factorization theorem for twist-two observables is also applicable at twist 
three. This leads to an expression that for intermediate qx agrees with the high-gr result 
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in its power behavior and in the form of the hard-scattering kernels, except for a term 
proportional to fi(x)D%(z). We find this partial agreement encouraging, but it does show 
that our candidate factorization formula at twist three is incorrect as it stands, and that 
a proper analysis will have to devote special attention to gluons with vanishing plus- or 
minus-momentum and to the precise form of soft factors. We emphasize that the correct 
description of F^^ h at low qx is a prerequisite for applying the method of Collins, Soper, 
and Sterman to resum large logarithms of Q 2 jq\- 

The structure function F^ T h ^ for a transversely polarized target presents a case 
where the low-q? calculation is of twist two, whereas the high-^ description is at the twist- 
three level. The explicit computations in [5-8] find exact agreement of the two descriptions 
at intermediate qT and thus validate both frameworks. One may expect that the same is 
true for F^ h+</,s) , which at low qr is described in terms of the Collins effect. 

Observables that are integrated over qj> are at times preferable to differential ones from 
an experimental point of view. We have shown that some of them have the added virtue of 
admitting a relatively simple description at the theory level, both for the complexity of the 
expressions and for the number of distribution and fragmentation functions on which they 
depend. With the power-counting behavior listed in table ^ one can readily determine to 
which region of qx a given integrated observable is primarily sensitive. The results for se- 
lected observables are given in table ||[ We find for instance that (( F^ 2 ^ h )) and (( F^ T ^ S )) 
receive leading contributions from both low and high qT- The integrated structure func- 
tion {{F^^ h )) is dominated by large qT, with contributions from the low-qy region being 
suppressed by M/Q. Conversely, both ((F^ h " 0s) )) and {{F^ hUs) » receive their dom- 
inant contributions from low qj>, whereas the high-^T domain is suppressed by M/Q. They 
are hence sensitive to the Sivers function in the first case, and to the transversity distri- 
bution and the Collins fragmentation function in the second. A suppression by M/Q may, 
however, not be sufficient to simply neglect the corresponding contributions in an analysis 
at experimentally achievable values of Q. 

A theoretically cleaner access to the high-gy region is through observables that are 
weighted with an appropriate power of qx/M. We find in particular that (( (qx/M) F™£ ^ h )) , 
(({q T /M) 2 F^ 2(t>h }}, and (({q T /M) 2 F^ s }} can be evaluated from the high-g T results 
alone, up to corrections of order M 2 /Q 2 , and that at the same accuracy one can extend 
the integration down to qT = 0. This leads to simple expressions, similar to the one for the 
integrated longitudinal structure function ((F UUL )). The observables {{{qr/M) )) 
and {{{qT / M) 2 F^ 2 ^ h )) are sensitive to the twist-two functions f\ and D\ and may for 
instance be useful for separating the contributions from different quark flavors, serving as 
complements to {{F UUT )). In contrast, {{{qT /M) 2 F^ T ^ S )) is sensitive to distribution and 
fragmentation functions of twist three. 

The weighted structure functions {{{q T /M)F^fp i ~ <l>s) )) and {{{q T /M)F s ^ 4 ' h+4 ' s) )) 
play a special role in this context. They receive leading-power contributions from both low 
and high qT and, as already pointed out in [12], lead to a deconvolution of the transverse- 
momentum integrals in the low-gr result. We argued that they should permit a description 
in terms of collinear functions of twist two and three, defined in the standard MS scheme. 
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In this description, the low-gy expression gives the Born- level result, whereas the high-gy 
calculation of [8,43,44] gives part of the a s corrections. If completed, such a description 
would provide a full NLO result in a s and be an extension to twist-three level of the stan- 
dard NLO calculation for {{Fjju.t)) within collinear factorization at twist- two accuracy. 
The leading-order expressions ( |7.9| ) and ( |7.10| ) for the weighted structure functions are 
analogs of the familiar tree-level formula = ^2 a xe a ff ( x ) (z). These expres- 

sions receive corrections from the high-gr region which are of leading power but suppressed 
by a s . 

Let us finally remark that the results we have discussed here carry over to the analo- 
gous observables in the Drell-Yan process and in e + e~ annihilation. The SIDIS structure 
functions F^^ h and F^ 2( t >h correspond for instance to the cos0 and cos 20 asymmetries 
in the angular distribution of the lepton pair in unpolarized Drell-Yan production, which 
have been measured [70] and given rise to several theoretical investigations, see e.g. the 
references in [9]. Furthermore, F^ 2 ^ h corresponds to a cos 20 asymmetry for two-pion 
production in e + e~ annihilation, which has been measured by BELLE [71] and provides 
the possibility for an independent determination of the Collins fragmentation function [72] . 
For a reliable extraction, our discussion of matching low- and high-transverse-momentum 
contributions should be of relevance. 
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A. Distribution functions at high px' Feynman versus axial gauge 



Our calculation in section 8.1 is done in axial gauge A ■ v = 0. Let us see how the same 
calculation proceeds in Feynman gauge. In this case one must explicitly take into account 
the gauge link IA in the definition fl5.3| ) of the correlation function $>(x,px), which consists 
of sections pointing along v and a transverse section at infinity. The detailed path of the 
gauge link reflects important physics, as shown for instance in [13,47-50,54]. 

Let us consider the correlation function $W(rc,p^) relevant for SIDIS, whose gauge 
link U + is closed at a~ = +oo. To evaluate the quark-to-quark contribution to the high- 
Pt behavior of §(x,pt) at leading order in 1/pr, one has to take into account the four 
diagrams shown in Fig. [?]. The graphs with eikonal lines are due to gluons coupling to the 
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(a) 



(b) 



l-p 








(< 






l-p 




(d) 



Figure 7: Diagrams for the calculation of the leading high-pr behavior of the quark-quark 
correlator $(x,pr) in Feynman gauge. 



gauge link in the operator tpj(0)U^ Q ^ V'i(C)- The corresponding Feynman rules read [24,55] 



ig t u v 



a, fi 



I ■ v + ie 



(A.l) 



where the sign of ie for the eikonal line corresponds to a gauge link pointing to a~ = +oo 
if one takes v ~ > 0. In cut diagrams one must take the conjugate of these expressions for 
vertices and propagators on the right of the final-state cut (indicated by the dashed lines 
in Fig. 0). 

After performing the integration over p~ using the 5 function in ( |8.3|) , we have for the 
respective diagrams 



^ [+]q (x,p T ) 
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*2 



x/ p z 

V 2 



(1-x) Vx/ [(/ -p) .« + i e ] [(/ -p) - ie] 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



with p and Z given in (|8.5[) and (8.6). Since p 2 is always spacelike according to (|8.7|) , we 
have omitted the ie in the quark propagators. 
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In the calculation using axial gauge only the first of the four diagrams contributes, but 
instead of g^ v in ( |A.2| ) we then have to take 



-*>«- = r - " ~ PTv - + " ~ PY v - + g -/'" ( ' - pT . 2 . (a.6) 

(I — p) ■ v [(I — p) ■ v] 2 

Each term in ( |A.6|) corresponds to one of the four diagrams in the Feynman gauge calcu- 
lation. The correspondence between the first term and diagram [|a in Feynman gauge is 
trivial. To establish the correspondence between the second term and diagram [?]b we use 
that 

**I(7-*)Js= Wji- m = -1*\, (A.7) 

where in the second step we have used that the twist-two part of the collinear quark correla- 
tor satisfies <3?2 / = ^2 ^+ ^ + = 0- In an analogous way one establishes the correspondence 
between the last two terms in ( A.6 ) with the respective contributions of diagrams 0c and 



0d in Feynman gauge. 

A mismatch between the expressions in ( |A.3| ) to ( |A.5| ) and the calculation in axial gauge 
is, however, the different treatment of the singularities at (l—p)-v = 0. The principal value 
prescription we employed when using the spacelike axial gauge of the original Collins- 
Soper paper [24] differs from the ie prescription for the different terms in the Feynman 
gauge calculation, which arises from the structure of the Wilson line IA + in the correlation 
function. We note that the integral in (A. 5) is actually not well defined as it stands, 



since the double pole at (I — p) ■ v = is pinched. The contribution from such unphysical 
poles must be absent in the physical cross section and should hence cancel between the 
distribution function, the fragmentation function, and the soft factor in the factorization 
formula ( |3.22 ). How to implement this by regulating the individual factors has so far not 



been addressed in the literature. It is also currently unknown if and how the principal- 
value prescription in axial gauge can be implemented in terms of Wilson lines for Feynman 
gauge. The discussion in [47] is for a light-cone rather than an axial gauge and hence does 
not contain the problematic term ( A.5|) . 



We remark that corresponding problems did not appear in the Feynman gauge calcu- 
lation of Ji et al. [26], where the vector v was chosen to be timelike. If we do the same in 
our context, then the calculations in Feynman and axial gauge exactly coincide. This is 
because (l—p)-v remains positive in (|A.3|) to ( A.5 ) according to ( 8.10|) , so that the singu- 



larity at (I— p)-v = is not reached in the loop integral. As a consequence, the particular 
regularization of the axial-gauge propagator does not influence our results of section || if 
we take v timelike. Likewise, there is then no contribution from transverse segments of the 
gauge link at infinity, which involve a 8 function in (I —p)-v. This is not implausible, since 
the distribution functions considered in section |8| are T-even and must in particular be the 
same for the gauge links U + and U~ . 

B. Integrated distribution functions and transverse-momentum cutoff 

In this appendix we derive Eq. ( |5.17| ) , which relates two different ways of regularizing the 
integral over the transverse-momentum-dependent distribution fi(x,p^). More precisely 
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we show that with \x = bo/b one has 

• f'OO 

d 2 p e* b - p f(p 2 ) = n / dp 2 J (bp) f(p 2 ) =77 I dp 2 f(p 2 ) + 0{b 2 \ 2 ) (B.l) 



for any function that can be expanded as 

f(p 2 )= C -i + C -i + C -i + ... (B.2) 

pZ p4 pD 

for p > A, where C2, C4, cq, etc. are constants. The intermediate scale A can be taken just 
large enough for (|B.2| ) to be valid, since corrections going for instance like M/A do not 
appear. The power corrections in ( |B. 1|) are understood as up to logarithms in b 2 X 2 . For 
ease of notation we have written p = \p\ and omitted the subscript T. 

To establish ( |B.1| ) we split the integrals into the regions p < A and p > A. In the first 
region we can write 

[ X dp 2 J (bp) f{p 2 ) = [ X dp 2 f{p 2 ) + 0(b 2 X 2 ) , (B.3) 
Jo Jo 

using that the Bessel function admits a Taylor expansion Jo(x) = 1 — \x 2 + . . . in even 
powers of x. In the region p > A we make use of the expansion ( p.2| ). Focusing first on the 
1/p 2 term, we write 

r°° dv 2 f°° dx f°° 

/ -Wo(M = 2/ — J Q (x) = -21n(6A) J (bX) + 2 / dxlnxJUx), (B.4) 
7a2 P 2 JbX x J bX 

where in the second step we have integrated by parts. We now use that dx lnx Ji(x) = 
In 2 — 7£ and Jq dx In x J\ (x) ~ y 2 for y — > 0, where the latter relation holds because 
J±(x) ~ x for x — > 0. Recalling that 60 = 2e~ 7B we obtain 

/7 |r J o(M = 2 N 2e ~ 7B ) - ln ( bA )] + o(& 2 a 2 ) = |r + °( b2x2 ) ■ (b.5) 

For the 1/p 4 term in ( |B.2 ) we use again integration by parts to write 

f°° ^ T (U \ 0* f°° dX T f \ J 0( &A ) h 2 [°°dx 
J\ 2 P JbX x A JbX x 

Since the integrand in the last term behaves like 1/x for x — > 0, we have 

d 4 Mb P ) = i + 0(6 2 A 2 ) = f d 4 + 0(b 2 X 2 ) . (B.7) 

A similar argument can be given for terms going like l/p 2n with n > 2, which completes 
the proof of Q. 

As an illustration of our result let us consider the simple form /(p 2 ) = l/(p 2 + A 2 ). 
The relevant integrals then are 

f°° dv 2 l'^ 2 dv 2 / li 2 \ 

1 1 J (bX) = 2K (bX), / _P =ln(l + ^). (B.8) 



/o p2 + A 2 uv v ur " 7 P 2 + A 2 V A 2 / ' 
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With the behavior Kq(x) = — Inx + lnbo + O^x 2 ) of the modified Bessel function at small x 
one readily finds that the relation ( |B.1| ) is satisfied. 

As we have seen in section 8.1, a logarithmic factor ln(£/p 2 ) appears in the explicit 
calculation for the high-pr behavior of distribution functions at order a s . One can easily 
repeat the above arguments for the case where the 1/p 2 term in (|B.2|) is multiplied by 



lnp 2 and the subleading terms by some power of Inp 2 . Using that dxihix) 2 J\{ 
(In 2 — 7_e) 2 one finds that (jB.ip holds without modification also in this case. 



C. One-loop expression of the soft factor 



In this appendix we show how to obtain the momentum-space expression flB.51 ) of the 



soft factor in the Collins-Soper factorization formula [24]. The corresponding expression in 
6-space is given in Eq. (7.22) of [24]. With our definition ( |3.7|) we obtain U{lj) from this 
by setting e = and omitting f c1 2 It e tb lT . The result is 

U{l T ) = -<iC F a s J -— p — (T (l;v) 

_ c F a * PV r di + u+v ~' v+ , =^', (CD 



^ 2 Jo ' [2{l+fv-/v+ +l 2 T ] 2 iy 

where on the last line we have indicated that for a spacelike gauge vector we need the 
principal value prescription to regulate the integral, given that v~ /v + < 0. In accordance 
with our footnote on page ||, the result of the integration is independent of v. Fourier 
transforming the result (C.l) to 6-space in 2 — e transverse dimensions, we obtain 



C F a s |£ f d 2 £ l T j lT . b 1 C F a s 



7T 2 J (27r)- £ 6 I 2 " 7T 



\nU 2 b 2 7re^)+- 



(C.2) 



in agreement with Eq. (7.23) in [24]. 
References 

[1] J. C. Collins, D. E. Soper, and G. Sterman, Transverse momentum distribution in Drell-Yan 
pair and W and Z boson production, Nucl. Phys. B250 (1985) 199. 

[2] R. N. Cahn, Azimuthal dependence in leptoproduction: a simple parton model calculation, 
Phys. Lett. B78 (1978) 269. 

[3] R. N. Cahn, Critique of parton model calculations of azimuthal dependence in 
leptoproduction, Phys. Rev. D40 (1989) 3107. 

[4] H. Georgi and H. D. Politzer, Clean tests of QCD in up scattering, Phys. Rev. Lett. 40 
(1978) 3. 

[5] X. Ji, J.-W. Qiu, W. Vogelsang, and F. Yuan, A unified picture for single transverse-spin 
asymmetries in hard processes, Phys. Rev. Lett. 97 (2006) 082002 friep-ph/0602239| . 



-67- 



[7; 



no 



[12 



[13 



[14 



[15 



[ie; 



X. Ji, J.-W. Qiu, W. Vogelsang, and F. Yuan, Single transverse-spin asymmetry in Drell-Yan 
production at large and moderate transverse momentum, Phys. Rev. D73 (2006) 094017 



| hep-ph/0604023 



X. Ji, J.-W. Qiu, W. Vogelsang, and F. Yuan, Single-transverse spin asymmetry in 
semi-inclusive deep inelastic scattering, Phys. Lett. B638 (2006) 178 | hep-ph/0604128 | 



Y. Koike, W. Vogelsang, and F. Yuan, On the relation between mechanisms for 
single-transverse- spin asymmetries, Phys. Lett. B659 (2008) 878 [{arXiv: 0711 . 0636 
[hep-phT| . 



D. Boer and W. Vogelsang, Drell-Yan lepton angular distribution at small transverse 



momentum, Phys. Rev. D74 (2006) 014004 [hep-ph/0604177| 



E. L. Berger, J.-W. Qiu, and R. A. Rodriguez-Pedraza, Transverse momentum dependence of 
the angular distribution of the Drell-Yan process, Phys. Rev. D76 (2007) 074006 



arXiv: 0708 .0578 [hep-ph] 



P. J. Mulders and R. D. Tangerman, The complete tree-level result up to order 1/Q for 



polarized deep-inelastic leptoproduction, Nucl. Phys. B461 (1996) 197 [tiep-ph/9510301]; 
Erratum-ibid. B484 (1997) 538. 

D. Boer and P. J. Mulders, Time-reversal odd distribution functions in leptoproduction, Phys. 



Rev. D57 (1998) 5780 [hep-ph/9711485 



D. Boer, P. J. Mulders, and F. Pijlman, Universality of T-odd effects in single spin and 



azimuthal asymmetries, Nucl. Phys. B667 (2003) 201 [hep-ph/0303034 



A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P. J. Mulders, and M. Schlegel, Semi-inclusive 
deep inelastic scattering at small transverse momentum, JHEP 02 (2007) 093 



[ hep-ph/0611265 



L. P. Gamberg, D. S. Hwang, A. Metz, and M. Schlegel, Light-cone divergence in twist-3 
correlation functions, Phys. Lett. B639 (2006) 508 |hep-ph/0604022| . 



EMC Collaboration, M. Arneodo et at, Measurement of hadron azimuthal distributions in 
deep inelastic muon proton scattering, Z. Phys. C34 (1987) 277; 

E665 Collaboration, M. R. Adams et al, Perturbative QCD effects observed in ^90 GeV deep 
inelastic muon scattering, Phys. Rev. D48 (1993) 5057; 

ZEUS Collaboration, M. Derrick et al., Inclusive charged particle distributions in deep 



inelastic scattering events at HERA, Z. Phys. C70 (1996) 1 [hep-ex/9511010]; 

HI Collaboration, C. Adloff et al., Measurement of charged particle transverse momentum 



spectra in deep inelastic scattering, Nucl. Phys. B485 (1997) 3 [iep-ex/9610006|; 
E665 Collaboration, M. R. Adams et al., Inclusive single-particle distributions and transverse 
momenta of forward produced charged hadrons in [ip scattering at 470 GeV, Z. Phys. C76 
(1997) 441; 

ZEUS Collaboration, J. Breitweg et al., Measurement of multiplicity and momentum spectra 
in the current and target regions of the Breit frame in deep inelastic scattering at HERA, 
Eur. Phys. J. Cll (1999) 251 |hep-ex/9903056| ; 

H. Mkrtchyan et al., Transverse momentum dependence of semi-inclusive pion production, 
Phys. Lett. B665 (2008) 20 [|arXiv : 0709 . 3020 [hep-ph] . 



[17] ZEUS Collaboration, J. Breitweg et al., Measurement of azimuthal asymmetries in deep 



inelastic scattering, Phys. Lett. B481 (2000) 199 [hep-ex/0003017 



- 68 - 



ZEUS Collaboration, S. Chekanov et at, Measurement of azimuthal asymmetries in neutral 
current deep inelastic scattering at HERA, Eur. Phys. J. C51 (2007) 289 | hep-ex/ 0608053 ] . 

[18] CLAS Collaboration, H. Avakian et al., Measurement of beam- spin asymmetries for deep 
inelastic 7r+ electroproduction, Phys. Rev. D69 (2004) 112004 |iep-ex/0301005| ]; 
HERMES Collaboration, A. Airapetian et al., Beam-spin asymmetries in the azimuthal 
distribution of pion electroproduction, Phys. Lett. B648 (2007) 164 |hep-ex/0612059]. 

[19] HERMES Collaboration, A. Airapetian et at, Observation of a single-spin azimuthal 
asymmetry in semi-inclusive pion electro-production, Phys. Rev. Lett. 84 (2000) 4047 
| liep-ex/991006^[ ; 

HERMES Collaboration, A. Airapetian et at, Single-spin azimuthal asymmetries in 
electroproduction of neutral pions in semi-inclusive deep-inelastic scattering, Phys. Rev. D64 
(2001) 097101 |hep-ex/0104005| ; 

HERMES Collaboration, A. Airapetian et at, Measurement of single- spin azimuthal 
asymmetries in semi-inclusive electroproduction of pions and kaons on a longitudinally 
polarised deuterium target, Phys. Lett. B562 (2003) 182 |hep-ex/0212039| ; 
HERMES Collaboration, A. Airapetian et al, Subleading-twist effects in single-spin 
asymmetries in semi-inclusive deep-inelastic scattering on a longitudinally polarized hydrogen 
target, Phys. Lett. B622 (2005) 14 |hep-ex/0505042( |. 

[20] HERMES Collaboration, A. Airapetian et at, Single-spin asymmetries in semi-inclusive 
deep-inelastic scattering on a transversely polarized hydrogen target, Phys. Rev. Lett. 94 
(2005) 012002 |hep-ex/0408013| ; 

HERMES Collaboration, M. Diefenthaler, HERMES measurements of Collins and Sivers 



asymmetries from a transversely polarised hydrogen target arXiv: 0706. 2242 [hep-ex] 



[21] COMPASS Collaboration, V. Y. Alexakhin et al, First measurement of the transverse spin 
asymmetries of the deuteron in semi-inclusive deep inelastic scattering, Phys. Rev. Lett. 94 
(2005) 202002 |hep-ex/0503002| ; 

COMPASS Collaboration, E. S. Ageev et al., A new measurement of the Collins and Sivers 
asymmetries on a transversely polarised deuteron target, Nucl. Phys. B765 (2007) 31 
[ hep-ex/0610068] ]; 

COMPASS Collaboration, A. Kotzinian, Beyond Collins and Sivers: further measurements of 
the target transverse spin- dependent azimuthal asymmetries in semi-inclusive DLS from 
COMPASS, larXiv : 0705 . 2402 [hep-ex] . 



[22] U. D'Alesio and F. Murgia, Azimuthal and single spin asymmetries in hard scattering 
processes, 



arXiv: 0712 .4328 [hep-ph] 



[23] A. Bacchetta, U. D'Alesio, M. Diehl, and C. A. Miller, Single-spin asymmetries: The Trento 
conventions, Phys. Rev. D70 (2004) 117504 ]hep-ph/0410050| |. 

[24] J. C. Collins and D. E. Soper, Back-to-back jets in QCD, Nucl. Phys. B193 (1981) 381. 

[25] J. C. Collins, Wha t exactly is a parton density?, Acta Phys. Polon. B34 (2003) 3103 
l hep-ph/0304122' l. 

[26] X. Ji, J. -P. Ma, and F. Yuan, QCD factorization for semi-inclusive deep-inelastic scattering 
at low transverse momentum, Phys. Rev. D71 (2005) 034005 Hhep-ph/0404183 [. 

[27] J. C. Collins, T. C. Rogers, and A. M. Stasto, Fully unintegrated parton correlation functions 
and factorization in lowest order hard scattering, Phys. Rev. D77 (2008) 085009 
[arXiv: 0708. 2833 [hep-ph] |. 



- 69 - 



P. M. Nadolsky, D. R. Stump, and C. P. Yuan, Semi-inclusive hadron production at HERA: 



The effect of QCD gluon resummation, Phys. Rev. D61 (2000) 014003 friep-ph/9906280| ; 
P. M. Nadolsky, D. R. Stump, and C. P. Yuan, Phenomenology of multiple parton radiation 



in semi-inclusive deep-inelastic scattering, Phys. Rev. D64 (2001) 114011 [tiep-ph/0012261| 



A. Kulesza, G. Sterman, and W. Vogelsang, Joint resummation in electroweak boson 



production, Phys. Rev. D66 (2002) 014011 fhep-pli/0202251 |. 

P. M. Nadolsky, Multiple parton radiation in hadroproduction at lepton hadron colliders, 



hep-ph/0108099 



S. Catani, D. de Florian, and M. Grazzini, Universality of non-leading logarithmic 
contributions in transverse momentum distributions, Nucl. Phys. B596 (2001) 299 



[hep-ph/0008184 



Y. Koike, J. Nagashima, and W. Vogelsang, Resummation for polarized semi-inclusive 
deep-inelastic scattering at small transverse momentum, Nucl. Phys. B744 (2006) 59 



[ hep-ph/0602188 ]. 

A. Weber, Soft gluon resummations for polarized Drell-Yan dimuon production, Nucl. Phys. 
B382 (1992) 63. 

P. M. Nadolsky and C. P. Yuan, Soft parton radiation in polarized vector boson production: 
Theoretical issues, Nucl. Phys. B666 (2003) 3 |hep-ph/0304001 1. 



J. C. Collins, Hard scattering in QCD with polarized beams, Nucl. Phys. B394 (1993) 169 



[ hep-ph/9207265 



X. Ji, J. -P. Ma, and F. Yuan, QCD factorization for spin- dependent cross sections in DIS 
and Drell-Yan processes at low transverse momentum, Phys. Lett. B597 (2004) 299 



[hep-ph/0405085l 



D. W. Sivers, Single spin production asymmetries from the hard scattering of pointlike 
constituents, Phys. Rev. D41 (1990) 83. 

J.-W. Qiu and G. Sterman, Single transverse spin asymmetries, Phys. Rev. Lett. 67 (1991) 
2264. 

A. Idilbi, X. Ji, J.-P. Ma, and F. Yuan, Collins-Soper equation for the energy evolution of 
transverse-momentum and spin dependent parton distributions, Phys. Rev. D70 (2004) 



074021 piep-ph/0406302 1. 

A. Mendez, QCD predictions for semi- inclusive and inclusive leptoproduction, Nucl. Phys. 
B145 (1978) 199. 

R. Meng, F. I. Olness, and D. E. Soper, Semi-inclusive deeply inelastic scattering at small qx, 



Phys. Rev. D54 (1996) 1919 |hep-ph/9511311 



M. Diehl and S. Sapeta, On the analysis of lepton scattering on longitudinally or transversely 
polarized protons, Eur. Phys. J. C41 (2005) 515 ftiep-ph/0503023 |. 



H. Eguchi, Y. Koike, and K. Tanaka, Single transverse spin asymmetry for large-px pion 
production in semi-inclusive deep inelastic scattering, Nucl. Phys. B752 (2006) 1 



| hep-ph/0604003 



H. Eguchi, Y. Koike, and K. Tanaka, Twist-3 formalism for single transverse spin asymmetry 
reexamined: Semi-inclusive deep inelastic scattering, Nucl. Phys. B763 (2007) 198 



[hep-ph/0610314 



- 70 - 



K. Hagiwara, K. Hikasa, and N. Kai, Time reversal odd asymmetry in semiinclusive 
leptoproduction in quantum chromodynamics, Phys. Rev. D27 (1983) 84. 



T. Gehrmann, Time-reversal- odd asymmetries at HERA, tiep-ph/9608469 



M. Ahmed and T. Gehrmann, Azimuthal asymmetries in hadronic final states at HERA, 



Phys. Lett. B465 (1999) 297 |iep-ph/9906503 |. 

A. V. Belitsky, X. Ji, and F. Yuan, Final state interactions and gauge invariant parton 



distributions, Nucl. Phys. B656 (2003) 165 friep-ph/020803S| 1. 

C. J. Bomhof, P. J. Mulders, and F. Pijlman, Gauge link structure in quark quark correlators 



in hard processes, Phys. Lett. B596 (2004) 277 |hep-ph/0406099| |. 

C. J. Bomhof, P. J. Mulders, and F. Pijlman, The construction of gauge-links in arbitrary 



hard processes, Eur. Phys. J. C47 (2006) 147 |iep-ph/0601171 |. 

J. C. Collins and A. Metz, Universality of soft and collinear factors in hard- scattering 



factorization, Phys. Rev. Lett. 93 (2004) 252001 ]hep-ph/0408249 |. 

F. Hautmann, Endpoint singularities in unintegrated parton distributions, Phys. Lett. B655 
(2007) 26 |hep-pti/0702196| ]. 

I. O. Cherednikov and N. G. Stefanis, Renormalization, Wilson lines, and 
transverse-momentum dependent parton distribution functions, Phys. Rev. D77 (2008) 



094001 |arXiv: 0710 .1955 [hep-ph][ ]. 

K. Goeke, A. Metz, and M. Schlegel, Parameterization of the quark-quark correlator of a 
spin- 1/2 hadron, Phys. Lett. B618 (2005) 90 |iep-ph/0504130| . 



J. C. Collins, Leading-twist single-transverse- spin asymmetries: Drell-Yan and deep-inelastic 
scattering, Phys. Lett. B536 (2002) 43 ftiep-ph/0204004 |. 



J. C. Collins and D. E. Soper, Parton distribution and decay functions, Nucl. Phys. B194 
(1982) 445. 



D. Boer, Theoretical aspects of spin physics, hep-ph/0312149 



A. A. Henncman, Scale dependence of correlations on the light-front, PhD thesis, Vrije 
Universiteit Amsterdam, 2005, 

http : //www.nikhef . nl/pub/ services/newbiblio/theses .php| . 

D. Boer, Sudakov suppression in azimuthal spin asymmetries, Nucl. Phys. B603 (2001) 195 



[ hep-ph/0102071 1. 

M. Anselmino, M. Boglione, U. DAlesio, A. Kotzinian, F. Murgia, and A. Prokudin, The role 
of Cahn and Sivers effects in deep inelastic scattering, Phys. Rev. D71 (2005) 074006 



[hep-ph/0501196l 



A. Mendez, A. Raychaudhuri, and V. J. Stenger, QCD Effects In Semiinclusive Neutrino 
Processes, Nucl. Phys. B148 (1979) 499. 

A. Konig and P. Kroll, A Realistic Calculation Of The Azimuthal Asymmetry In 
Semiinclusive Deep Inelastic Scattering, Z. Phys. C16 (1982) 89. 

J. Chay, S. D. Ellis, and W. J. Stirling, Azimuthal asymmetry in lepton-photon scattering at 
high energies, Phys. Rev. D45 (1992) 46. 

K. A. Oganessyan, H. R. Avakian, N. Bianchi, and P. Di Nezza, Investigations of azimuthal 
asymmetry in semi-inclusive leptoproduction, Eur. Phys. J. C5 (1998) 681 [ hep-ph/970934"^ |. 



- 71 - 



[64] M. Anselmino, M. Boglione, A. Prokudin, and C. Turk, Semi-inclusive deep inelastic 

scattering processes from small to large Py, Eur. Phys. J. A31 (2007) 373 | hep-ph/ 060628^] . 

[65] A. M. Kotzinian and P. J. Mulders, Longitudinal quark polarization in transversely polarized 
nucleons, Phys. Rev. D54 (1996) 1229 |iep-ph/951142C| ]. 

[66] J. Chay and S. M. Kim, Azimuthal correlation in lepton hadron scattering via charged 
weak-current processes, Phys. Rev. D57 (1998) 224 |hep-pti/9705284 ]. 

[67] J. C. Collins, Fragmentation of transversely polarized quarks probed in transverse momentum 
distributions, Nucl. Phys. B396 (1993) 161 frep-ph/9208213; ]. 

[68] S. Meissner, A. Metz, and K. Goeke, Relations between generalized and transverse momentum 
dependent parton distributions, Phys. Rev. D76 (2007) 034002 |iep-ph/0703176|1 . 

[69] X. Artru and M. Mekhfi, Transversely polarized parton densities, their evolution and their 
measurement, Z. Phys. C45 (1990) 669. 

[70] NA10 Collaboration, S. Falciano et al., Angular distributions of muon pairs produced by 194 
GeV/c negative pions, Z. Phys. C31 (1986) 513; 

NA10 Collaboration, M. Guanziroli et al., Angular distributions of muon pairs produced by 
negative pions on deuterium and tungsten, Z. Phys. C37 (1988) 545; 
J. S. Conway et at, Experimental study of muon pairs produced by 252- GeV pions on 
tungsten, Phys. Rev. D39 (1989) 92; 

FNAL-E866/NuSea Collaboration, L. Y. Zhu et al., Measurement of angular distributions of 
Drell-Yan dimuons in p + d interaction at 800 GeV/c, Phys. Rev. Lett. 99 (2007) 082301 
| liep-ex/0609005 l. 

[71] BELLE Collaboration, K. Abe et al., Measurement of azimuthal asymmetries in inclusive 
production of hadron pairs in e + e~ annihilation at BELLE, Phys. Rev. Lett. 96 (2006) 
232002 |hep-ex/0507063 ]. 

[72] D. Boer, R. Jakob, and P. J. Mulders, Asymmetries in polarized hadron production in e + e~ 
annihilation up to order 1/Q, Nucl. Phys. B504 (1997) 345 |hep-ph/970228l| ]. 

[73] D. Binosi and L. Thcussl, JaxoDraw: A graphical user interface for drawing Feynman 
diagrams, Comput. Phys. Commun. 161 (2004) 76 [ tiep-ph/0309015 |. 



- 72 - 



